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Abstract 

In this article, we study parameter-dependent systems of reaction-diffusion equations on the 
sphere, which are equivariant under the group 50(3) of all rigid rotations on a sphere. Two 
main types of spatial-temporal patterns that can appear in such systems are rotating waves 
(equilibrium in a co-rotating frame) and modulated rotating waves (periodic solution in a co- 
rotating frame). The transition from rotating waves to modulated rotating waves on spherical 
domains is explained via a supercritical Hopf bifurcation from a rotating wave, 5'0(3)-symmetry 
and finite-dimensional equivariant center manifold reduction. The Baker-Campbell-Haussdorff 
formula in the Lie algebra so(3) is used to get reduced differential equations on so(3), a formula 
for a primary frequency vector, as well as a formula for the periodic part associated to any 
modulated rotating wave obtained by a supercritical Hopf bifurcation from a rotating wave. As 
a consequence, there are three types of motions for the tips of these modulated rotating waves. 
In the resonant case (with two parameters) , we obtain that the primary frequency vectors of a 
branch of these modulated rotating waves are generically orthogonal to the frequency vector of 
the initial rotating wave undergoing Hopf bifurcation. 

Keywords : equivariant center manifold, rotating wave, modulated rotating wave, Hopf bifurca- 
tion 

1 Introduction 

The main motivation of this article is the presence of spiral waves in excitable media, especially in 
cardiac tissue. Spiral waves arise as stable spatio-temporal patterns in various chemical, physical 
systems and biological systems, as well as numerical simulations of reaction-diffusion systems on 
excitable media with various geometries. Excitable media are extended non-equilibrium systems 
having a uniform rest state that is linearly stable but susceptible to finite perturbations. Spi- 
ral waves have been observed experimentally, for instance, in catalysis of platinum surfaces j35j . 
Belousov-Zhabotinsky chemical reactions |25| I29j . Rayleigh-Bernard convection . slime-mould 
cells |13| and the most important, cardiac tissue [S]. Numerical simulations of the spiral waves have 
been done by 121 , for example. 
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It is now believed that spiral and scroll waves that appear in the heart muscle can lead to cardiac 
arrhythmias (abnormal rhythms in the heart), giving rise to atrial fluttering or ventricular fibril- 
lation. In normal hearts cardiac arrhythmias are rare, but in diseased hearts cardiac arrhythmias 
can become more common. For example, if chambers of the heart become abnormally large, they 
are susceptible to serious arrhythmias in which waves are believed to circulate in a fashion that is 
similar to the circulation of the Belousov-Zhabotinsky waves in a chemical medium. Real human 
hearts are enormously complex three-dimensional structures. In this article, we assume that the 
geometry of the excitable media is a sphere, which in the case of cardiac tissue is clearly an ap- 
proximation. 

In the planar case, a rigidly rotating spiral wave is an example of wave pattern rotating around a 
center and being well approximated by an Archimedean spiral wave far from rotation center. Near 
the rotation center, there is a core region of the spiral wave, where the front of the wave has a tip, 
whose structure is considered to be the most important in understanding the behavior of the whole 
spiral wave |261 127j . Barkley ^ was the first who performed a numerical linear stability analysis 
for the basic-time periodic spiral wave solution in a reaction-diffusion system on the unbounded 
plane and showed evidence of a Hopf bifurcation. In particular, a simple pair of eigenvalues was 
shown to cross the imaginary axis while three neutral eigenvalues lie on the imaginary axis and 
the remainder of the spectrum is bounded into the left-half plane. Later, using an ad hoc model, 
Barkley [2] was the first to realize the key importance of the group SE{2) of all planar translations 
and rotations in describing the dynamics and bifurcations of planar spiral waves. 
It is well known now that the tip of the spiral wave rotates steadily or meanders or linearly drift 
in plane jl4| From a mathematical point of view, rigidly rotating spiral waves are examples 
of rotating waves, meandering spiral waves are examples of modulated rotating waves and linearly 
drifting spiral waves are examples of modulated rotating waves. 

The first rigorous mathematical theory of the planar spiral waves was done by Wulff 52 . In her 
thesis, WulfF studied the external periodic forcing of rotating waves which leads to modulated ro- 
tating waves or modulated travelling waves, and also proved an S'£'(2)-equivariant Hopf theorem 
for the bifurcation from rotating waves to modulated (rotating or travelling) waves in autonomous 
systems. The external periodic forcing of rotating waves was studied using a contraction mapping 
theorem on scales of Banach spaces, and the proof for the S'i?(2)-equi variant Hopf theorem for 
rotating waves was based on Liapunov-Schmidt reduction on scales of Banach spaces. In both 
cases, it is shown that modulated travelling waves emanate if the rotation frequency is a multiple 
of the external frequency, respectively the modulus of the eigenvalue leading to Hopf bifurcation. 
The main difficulty comes from the fact that the group SE{2) is noncompact and the action of 
SE{2) on the usual spaces of functions is not smooth. The proofs in [32 are based on the basic 
assumption that the linearization at the rotating wave in the co-rotating frame does not exhibit 
continuous spectrum near the imaginary axis. 

Later, Sandstede, Scheel and Wulff [321 proved a finite-dimensional center bundle reduction theo- 
rem near a relative equilibrium Guq of an infinite-dimensional vector field on a Banach space X 
on which acts a finite-dimensional Lie group (not necessarily compact). This generalizes Krupa's 
results on bifurcation from relative equilibria (22] , from compact groups to noncompact groups and 
from finite dimensions to infinite dimensions. Using the results of jl^]) the Hopf bifurcation from 
one-armed rotating spiral wave to meandering waves can be studied jl4| I45j . 

In the case of a one-armed rotating spiral wave, the relative equilibrium SE{2)uq is diffeomorphic to 
SE{2), which is diffeomorphic to C x S^. Near Hopf bifurcation, the reduced differential equations 
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on the center bundle SE{2) x C are given by 



V = 



^ = 

q = 



(1.1) 



where F''^(0,0) = ujrot, F'^{0,0) = 0, /(0,0) = 0. The rotating wave uq corresponds to (; = at 
A = and to the solution (0,uJrott) for the first two differential equations in (|l.lj) . 
In case of a supercritical Hopf bifurcation we have dqF'^(0, 0) = 11^2, Re {dgF^{0, 0)) 7^ 0. It follows 
that for any small A > there is a meandering spiral wave which becomes a drifting linear wave 
if wi(A) = A;a;2(A), for some k £ Z, where u;2(A) is the frequency that appears due to the Hopf 

bifurcation and a;i(A) = ^"^^^^^^ Jo^^^^ F'^{q{t, A), A) dt, and q{t, A) is the |^^^^-^| periodic solution of 
the third differential equation in (|1.1() that appears due to the supercritical Hopf bifurcation for 
g = at A = 0. Similar results for m-armed spiral waves were also obtained. 

In |53j . a G-equivariant semilinear system of parabolic equations (where G is a finite-dimensional 
possibly non-compact Lie group) is studied. In particular, the periodic forcing of relative equilib- 
ria and resonant periodic forcing of relative equilibria to relative periodic orbits, as well as Hopf 
bifurcation from relative equilibria to relative periodic orbits are treated using Liapunov-Schmidt 
reduction. Resonant drift phenomena are also studied. Then, these results are applied to the planar 
spiral waves. 

For the skew-product finite-dimensional system of differential equations on the center manifold near 
a relative equilibrium , the normal form method which further simplifies the system is presented 
in |inj . The normal form for the case G = SE{N) is obtained. Then, the known results regarding 
the meandering and drifting of planar spiral waves are recovered, as well as new results regarding 
the relative homoclinic and heteroclinic trajectories to relative equilibria of S'£'(2)-actions. 
In Scheel shows that for a large class of reaction-diffusion systems on the plane, m-armed spi- 
ral waves bifurcate from a homogeneous equilibrium when the latter undergoes a Hopf bifurcation. 
This was done using spatial dynamics. 

All previous results are valid for planar spiral waves. The interest to consider spiral waves on non- 
planar surfaces is motivated by the applicability to problems in physiology, biology and chemistry. 
Therefore, the study of spiral waves by experiments and numerical simulations of reaction-diffusion 
systems on a sphere have recently been undertaken. In the case of spiral waves on a sphere, the 
dynamics is expected to be quite different because any spiral wave starting from a rotating center 
cannot end at a point. The number of tips of a wave front cannot be odd, and therefore, the 
dynamics of spiral waves may acquire a new feature qualitatively different from the planar case. 
The dynamics of spiral waves in an excitable reaction-diffusion systems on a sphere was numeri- 
cally investigated by j^HSHlEni) and [3^ who employ a spectral method using spherical harmonics 
as basis functions. Maselko I32j. as well as Maselko and Showalter jHSl performed experiments 
with Belousov-Zhabotinsky chemical waves propagating on the surface of a sphere. They observed 
that a spiral wave winds outward from a meandering source at the north pole and undergoes self- 
annihilation as it winds into itself at the south pole. 

In the evolution of spiral waves on a circular domain and on a spherical surface is studied by 
numerical integration of a reaction-diffusion system. Two different asymptotic regimes are observed 
for both domains. The first regime is a rigid rotation of an excitation wave around the symmetry 
axis of the domain. The second one is a compound rotation including a drift of the rotation center 
of the spiral wave either along the boundary of the disk or along the equator of the sphere. In this 
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case the shape of the wave and its rotation velocity are periodically changing in time. Simplified 
analytical estimates are presented to describe the rigid rotation. 

In jl8j . numerical integration of an excitable reaction-diffusion system on a sphere is presented. The 
evolution of counter-rotating double spiral waves on the sphere is studied and it is shown that the 
tips of the spiral can either perform a meandering motion or rigidly rotate around a fixed center, 
depending on the system control parameter. It is observed that the rotation of the spiral wave on 
spherical surface is similar to that obtain on the planar surface, except that in the absence of the 
boundary on a spherical surface some parts of the wave can undergo self-annihilation in contrast 
to the spiral wave behavior on bounded planar surfaces. 

In |56j , the evolution of spiral waves on a spherical surface is studied by integration of a reaction- 
diffusion system with a global feedback. It is shown that depending on intensity, sign, and/or 
time delay in the feedback loop a global coupling can be effectively used either to stabilize the rigid 
motion of a spiral wave or to completely destroy spiral waves and to suppress self-sustained activity 
in a confined domain of an excitable medium. 

In ) the dynamics of chemical spiral waves in an excitable reaction-diffusion system on a sphere is 
numerically investigated employing a spectral method using spherical harmonics as basis functions. 
Different types of spiral waves -symmetric or antisymmetric source-source or nearly antisymmet- 
ric source-nonsource— have been obtained depending on whether the medium is homogeneous or 
inhomogeneous, and it has been observed that the tips can either rotate steadily or change their 
shapes. 

The influence of the topological constraints and the inhomogeneity in the excitability on the ge- 
ometry and dynamics of spiral waves on a thin spherical shell of excitable media are presented in 
[3]. Also, rigidly rotating waves on spherical domains have been studied using kinematical theory 

by [in,2a..34 . 

In , the geometrical stability of the symmetric counter-rotating spiral waves propagating on the 
unit sphere is studied. By the use of the eikonal equations, it is demonstrated that these solutions 
are stable under small perturbations normal to the wave front lying on a unit sphere. 
In [20], the authors showed that stationary rotating solutions on a sphere of the eikonal equation 
under some boundary conditions on tips must be symmetric with respect to the equator with spiral 
winding out from source tips at polar points. 

In j34j . using the eikonal approximation to a reaction-diffusion system on a sphere, the authors 
prove existence of a class of counter-rotating double spiral solutions which are highly asymmetric 
with respect to the equator. They also derive a power law, linking the angular rotation of the spiral 
waves with the velocity of plane waves in medium. 

In Renardy considered bifurcations from rotating waves of semi-linear equations that are 
equivariant under a general compact Lie group and applied his results to the Laser equations. His 
results do not cover the resonance case, which will be covered in our paper. The theorems were 
proved using a generalized implicit function theorem on scales of Banach spaces. 
In ^4Qj, Rand examined modulated rotating waves in rotating fluids and applied his results to the 
Taylor-Couette problem, where modulated rotating waves, so called modulated wavy vortices oc- 
cur. 

In this article we have studied only the trivial isotropy case, because as far as we are aware, there 
have been observed no m-spiral waves (m > 1) on spherical surfaces. Also, in [S] it was numerically 
verified that there is a critical size of the sphere below which self-sustained spiral waves cannot 
exist. Therefore, we use a sphere of an arbitrary, but fixed radius r. 
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In this article, we study non-resonant and resonant supercritical Hopf bifurcation from a rotating 
wave. The relation between spiral waves and rotating waves, respectively modulated rotating waves 
will be presented somewhere else. While we have concentrated on the case of Hopf bifurcation, it 
should be mentioned that we would obtain similar results for periodic forcing of a rotating wave. 
In Section [21 we define the representation of the Lie group 5*0 (3) on the usual function space. It is 
a smooth unitary representation on the fractional spaces relative to the appropriate sectorial oper- 
ator. We also associate an S'0(3)-equivariant semiflow to the reaction-diffusion system on a sphere, 
and we recall the notion of tip position function. In Section El we recall the notions of relative 
equilibrium, relative periodic orbit, rotating wave and modulated rotating wave. We also define the 
notions of frequency vector of a rotating wave and primary frequency vector of a modulated rotating 
wave, and present some of their properties. Section deals with a supercritical Hopf bifurcation 
from a rotating wave. Decomposition of a modulated rotating wave obtained by Hopf bifurcation 
into a primary frequency part and the associated periodic part is proved. In Sectional the form of 
the reduced differential equations on so(3), as well as formulae for a primary frequency vectors and 
the associated periodic part are obtained using Baker-Campbell-Hausdorff formula in the Lie alge- 
bra so{3), the properties of the exponential map e: so(3) SO{3) and the properties of the adjoint 
representations of 5*0(3) and so{3). As a consequence, there are three types of motions for the tips 
of the modulated rotating waves obtained by a Hopf bifurcation from a rotating wave. Namely, the 
tip can either quasi-periodically meander on the sphere, such that the primary frequency vector is 
of order 0(1) near the frequency vector of the initial rotating wave, or quasi-periodically meander 
very slowly (order 0(\/A)) about the primary frequency vector, which can be orthogonal or not to 
the frequency vector of the initial rotating wave undergoing Hopf bifurcation. In the orthogonal 
case, the tip motion can be approximated with a very slow drift (order 0(\/A)) along an equator 
of the sphere. Section El presents some examples that illustrate the theoretical results obtained 
in Sections jll and O In Section [71 the resonant Hopf bifurcation from a rotating wave with two 
parameters is studied. Near the resonant Hopf bifurcation, we proved that there exists generically 
a unique branch of modulated rotating waves with primary frequency vectors orthogonal to the 
frequency vector of the rotating wave undergoing Hopf bifurcation. Then, Section |H1 deals with 
some numerical results that illustrate the theoretical results obtained in Sections 01 and O The 
proofs of all theorems are presented in Sectional Three appendices are included, presenting BCH 
formula in so(3), some properties of the exponential map e: so(3) — > 50(3) and of the adjoint 
representations of 50(3) and so(3), the equivariant center manifold reduction theorem and some 
computations done in Section |H1 



2 Reaction-diffusion Systems on Sphere rS^ 

Let r > and be the unit sphere in M'^. We consider a reaction-diffusion system of the form 





U = (lii, U2, ■■■ 

I di ... 
D = ■■ ■. 




Uiv) : M X rS^ with N > 1, 

M 

'■ with di > for i = 1, 2, . . . , N the diffusion coefficients, is Laplace- 
dN J 
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Beltrami operator on rS^ and F = (Fi, Fa, . . . , Fat) : . Also, DAsu 



We study the reaction-diffusion system t\2.1^ on the function space 



/ diAsUi \ 

d2AsU2 

\ (InAsun J 




if > for i = 1, 2, . . . , N; 
if there exists i £ {1, 2 ... , N} such that di = 0. 

Let a E (i, 1) and Y° = D{{I - DAs)°') be the fractional spaces of Y relative to the sectorial 
operator —DAs- Using ^1 1^ OHl 113 IS2] j to the reaction-diffusion system ()2.1() we associate 
a local semiflow on an open set of Y. Let F: ^ be a C*^+2 function such that F(0) = 0, 
where < A; < oo. Then the reaction-diffusion system H2.1|) defines a sufficiently smooth local 
semiflow $ on the function space Y°. The flow is (7^+^ if y = L2(rS2,M^) and CMf Y = 
H2(rS2,M^). Namely [SHI, for any uq G Y", let u{t,uo) be the sufficiently smooth solution {C''+^ 
or C'^ respectively) of the initial value problem given by the reaction-diffusion system (|2.1|) and by 
the initial condition u(0) = uq, defined on the maximal interval of existence I{uq) = [0,to('"o))- 
Let W = {{t,uo) e [0,oo) X Y" I t E /(no)}. Then, the local semiflow 1^ ^ Y" is defined 
by $(t,uo) = u{t,uo), for any (t, mq) G W. Any reaction-diffusion system of the above type has 
the property that is S'0(3)-equivariant, where the precise meaning of 50(3)-equivariance is defined 
below. 

Definition 2.1. The representation T of SO {3) on Y is the function T : S0{3) GL(Y), defined 
by 

T{A)u{x) = u{A'^x) where A £ S0{3),u G Y, x G rS^. (2.2) 

The linear action 9: S0{3) x Y ^ Y associated to the unitary representation T ( |22[ I23[ 130 1 
I1I1E2]) is defined by 

9{A, u){x) = u{A-^x) where A G S'0(3), n € Y and x € rS^. (2.3) 

We denote 9{A, u) = Au. From now on, when we talk about 5'0(3)-equivariance, we mean equivari- 
ance with respect to the action 9 defined in 1)2. 3|) . Using [Tj |^ 0H1 IMl EOl > we obtain the following 
two propositions: 

Proposition 2.2. ^ The restriction of the representation T to Y", where a G (^,1); is smooth. 

Proposition 2.3. \^ Let F:R^ ^ be a C^^"^ function such that F(0) = and 1 < A; < oo. 

The local semiflow $ is SO{3)-equivariant with respect to the action 9 restricted to Y" 



For any no € Y° we follow the time evolution of uq under the local semiflow $ using the concept 
of tip motion of uq. 

Definition 2.4. For < A; < oo, a function xup'. Y" — > rS^ that is SO{3)-equivariant is 
called a tip position function. By the tip of u £ Y", we understand the point xtip{u) S rS^. Also, 
for any uq E Y", xtip{^{t, no)) with t G [0, cxd) is called the tip motion of uq. If the function xup is 
defined only on an open subset ofY'"', then the function xup is called a local tip position function. 

The tip position function is usually used to follow the time evolution of no under the local 
semiflow when $(t, no) is a rotating wave or a modulated rotating wave. 
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3 Rotating Waves and Modulated Rotating Waves on rS^ 



We consider that north hemisphere of is the set = {(x, y, z) G | 2 > or z = 0, x G 
[—1, 1), y G [0, 1]}, and south hemisphere is the set S = S"^ \ N. 

For any G-equivariant dynamical system, where G is a Lie group, we can talk about the concepts 
of relative equilibrium and relative periodic orbit. If SO{n) C G for some integer n > 2, we can 
also talk about rotating waves and modulated rotating waves. 

Definition 3.1. fTTl US E^/ Let uq G Y" be such that the stabilizer of uq is — I^. The orbit 
group of Uq is called a relative equilibrium for \2.1\) if there exists a matrix Xq G so(3) such that 

Uq) = e^°*Uo for any t > 0. (3.1) 

Sometimes, when S0{3)uo is a relative equilibrium, we call uq a relative equilibrium. If Xq ^ O3, 
then any solution of the reaction- diffusion system \2. 1\) of the form A^{.,uq), where A G S'0(3) is 
called a rotating wave for \2.1\) . 

Because the action 9 restricted to Y" is smooth, then the relative equilibrium from Definition 
13. II is a sufficiently smooth manifold in Y" diffeomorphic to S'0(3) since S„q = I^. 

Definition 3.2. Let ^{.,Auq) be a rotating wave as in Definition Vj.lV where A G S'0(3). Then, 
> > 

the vector AXqA~^ is called the frequency vector of the rotating wave. If the vector jy^AXqA~^ 
is in the north hemisphere, then loq = \Xq\ is called the frequency of the rotating wave. Otherwise, 
—ujQ is called the frequency of the rotating wave. 

Then, the rotating waves in SO{3)uq have their anti-symmetric matrices associated with the 
frequency vectors on the same adjoint orbit. The following property shows the importance of the 
frequency vector associated to any rotating wave: 

Proposition 3.3. ^ For a rotating wave mq); the tip motion xtip{^{t,uo)) is a circle on the 
sphere rS^ with the center on the line having the direction of the frequency vector of ^{t,Uo), and 
this is independent of the choice of the tip position function. 

Definition 3.4. lTli\EEi\SBl Let uq G Y" be such that the stabilizer of uq is Jjjio — -^3- The set 

defined by mq) | A G 5*0(3), t G [0,oo)} is called a relative periodic orbit for j() if it is 

not a relative equilibrium and there exist a number T > and a matrix Xq G so(3) such that 

$(T,uo) = e^^^uo and ^{t,uo) ^ 50(3)uo for any t G (0,r). (3.2) 

If \Xq\ T 7^ 2A;7r for any G Z, then any solution of the reaction- diffusion system h2. 1)) of the form 
A^{.,uq), with A G 50(3) is called a modulated rotating wave for \2.1]) . 

If |Xo| = ^ for some k £ Z, then $(t, mq) is a T-periodic solution of the reaction-diffusion 
1)2.11) . Because the action of SO (3) on Y*^ is smooth, then the relative periodic orbit defined 
in Definition 13.41 is a sufficiently smooth manifold in Y". There are two properties related to 
modulated rotating waves, namely: 

Remark 3.5. f4^I 

1. *(ri,uo) = e^o^iiio if and only if Ti G TZ. 
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2. There are periodic solutions of the reaction- diffusion system \2.1\) that are modulated rotat- 
ing waves in the sense of Definition Then their period is an integer multiple of the 
corresponding T from Definition \3.4\ and \Xq\ T S 27rQ. 

Definition 3.6. Let $(., Auq) be a modulated rotating wave as in Definition \°j.4\ where A £ S0{3). 

Then, the vector AXqA~^ is called a primary frequency vector of the modulated rotating wave and 

the positive number T is called the relative period of the modulated rotating wave. If the vector 
> 

jY^AXqA"^ is in the north hemisphere, then ojq = \Xo\ is called the primary frequency of the 
modulated rotating wave. Otherwise, —loq is called the primary frequency of the modulated rotating 
wave. 

A primary frequency vector Xi of a modulated rotating wave ^{t, Auq) is unique up to an integer 
multiple of j^^Xi. Then the modulated rotating waves in 5'0(3)$(., uq) have the anti-symmetric 
matrices associated with primary frequency vectors on the same adjoint orbit (if we consider primary 
frequency vectors to be unique up to a multiple of /c G Z of the corresponding unit primary 
frequency vector). If ^mq) is a T-periodic solution of the reaction-diffusion system such 
that ^{T,Auq) = Auq and ^{t,AuQ) ^ S0{3)uq for all t G (0,T), then a primary frequency vector 
of ^{t,Auo) can be any vector Xi G with \Xi\ = for A; G Z. The following property shows 
the importance of the primary frequency vector associated to any modulated rotating wave: 

Proposition 3.7. For a modulated rotating wave uq), the tip motion xtip{^{t,UQ)) has the 
property that xtip{^{kT,uo)) for k £ Z, are points of a circle on the sphere rS^ with the center on 
the line having the direction of the primary frequency vector of ^{t, uq); this is independent of the 
choice of the tip position function. 



4 Hopf Bifurcation from Rotating Waves to Modulated Rotating 
Waves on rS^ 

In this section, we will consider a supercritical Hopf bifurcation from an equilibrium. Because the 
bifurcating periodic solution has amplitude of order \/A, the following definition will simplify the 
language. 

Definition 4.1. Let M be a smooth manifold, X a normed space or the empty set, p > an 
integer and Y : X x [0, Aq) x W ^ M for Aq > small. We say that the function Y is CS on 
X X [0, Ao) X W if the function Z: X x[{), eo) xW ^ M defined by 

Z{x,e,fi) = y(a;,e^/i), 

is smooth on X x [0, eo) x , where eq = \/Ao- 

We say that Y is C^-CS if Z is , where k e Z, k > 1. We say that Y is sufficiently CS if Z is 
sufficiently smooth. 

Let us consider the following reaction-diffusion system 

— {t, x) = DAsu{t, x) + F{u{t, x), A) on rS^, (4.1) 
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( di ... 



d,- > 0, i = 1, 2, 



u = (ui, n2, . . . , uat) : M X rS^ ^ with N > I, D = 

\ ... djv 

. . . , TV the diffusion coefficients and F = (Fi, F2, . . . , Fat) : x M ^ a C''+2 function with 
3 < A; < CO. We study the reaction-diffusion system (|4.1j) on the function space Y defined in 
Section [21 Let mq £ be a relative equihbrium that is not an equilibrium for (|4.ip at A = and 
such that the stabilizer of uq is = J3. Let mqiO) = e"'''''*no. Consider L the linearization of 
the right-hand side of (|4.1|) with respect to the rotating wave ?io,0) = e-^°^uo at A = in the 
co-rotating frame, that is 

L = DAs + DuF{uo,0) - Xo. 

Suppose that: 

1. (j{L)n{z £ C I Re (z) > 0} is a spectral set with spectral projection P*, and dim(i?(P*)) < 00; 

2. the semigroup e^* satisfies |e^*|_R(i-p,) | < Ce~^°* for some /3o > and C > 0. 

Theorems IB. II (see Appendix E| can be applied. Then, there exist sufficiently smooth functions 
Xg : K X K — > so(3) and Xj\f : K x 1^ — > K such that any solution of 

A = AXG{q,X), /4 2^ 

q = XN{q,X), ^ ■ ^ 

on 50(3) X corresponds to a solution of the reaction-diffusion system (|4.H) on Af^"(A) under the 
diffeomorphic identification for |A| small. Also, Xg{0,0) = Xq, X]\f{0,0) = and a{DuX]\f{0,0)) = 
a{Q^:L\v^), where Q=k is the projection onto K along Tuq{SO{3)uo). An example of tip position 
function is xup'. M™(A) x M — > rS^ defined by xtip{A'i^{q), X) = Axq, where xq € rS^ is fixed. 
Also, we have the following result. Suppose we omit the parameter A. Let $(.,ui) be a modulated 
rotating wave as defined in Definition 13.41 with ui E M^" and = I^, where M^" is defined in 
Theorem IB. II Suppose that ^{.,uq) corresponds to and $(.,mi) corresponds to qi{t), where 
and qi{t) are solutions of the second differential equation in (|4.2j) . Using 9. .1"7| BHj. we have: 

Proposition 4.2. ^ 

1. $(.,Mo) is orhitally stable (respectively unstable) if is stable (respectively unstable) in the 
second differential equation of i4-^ - 

2. $(.,ni) is orbitally stable (respectively unstable) if qi{t) is stable (respectively unstable) in the 
second differential equation of i4-^ - 

Let dim(i?(Fi<)) = 5. Suppose that a supercritical Hopf bifurcation with eigenvalues ^iutif 
takes place in the second differential equation of 1)4. 2() in K at g = for A = |22 . Namely, 

L Xiv(0,0) =0; 

2. DqXj\f{0,0) has eigenvalues ^iujuf'-, without loss of generality, we assume that XAr(0, A) = 
for |A| small; 

3. -DgXAr(0, A) has the eigenvalues a(A)ibi(u;fej/+/3(A)) with a(0) = /?(0) = such that a' (0) > 0; 
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4. the branch of periodic solutions q{t,X) bifurcating from q = satisfies q{t,\) = 0(\/A) for 
A > small; 

Let q{t, 0) = for all t G [0, cxd). Let T(A) = |^ be the period of the solution q{t, A) near q = 
that appears for A > small due to the supercritical Hopf bifurcation, where ujx = Ubif + 0(A) for 
A > small. Let us define 



iXG{q{t,X)A) if A > and t G [0, cx)) 



[Xo if A = and t G [0,oo), 



where XG{q,X) is defined in Theorem IB.lf see Appendix IB)) . 

Since q{t,X) = \/l^r{t,X) is sufficiently CS and XG{q,X) is sufficiently smooth, it follows that the 
function X'-' is sufficiently CS. 
Writing q{t, X) = VXr{t, A), we have 



X^{t,X) = Xo{t,X)X^ + ^/Xxl{t,X)Xl+^X2it,X)X2, 

xo{t,X) = |Xo| + \/Axoi(t) + Aa;o2(t, A) 



(4.4) 



for t G [0, oo) and A > small. The functions X'^{t,X), xo{t,X), xi(t. A), X2{t,X) are |^-periodic 
in t for A > small. 

Let A{t, A) for t G [0, cxd) and A > small be the solution of the initial value problem 

A = AX^itA), 

A{0) = h, ^^-^^ 

where X^{t, A) is defined in (|4.3)) . 

Lemma 4.3 (Decomposition of A{t,X)). Suppose that the assumptions made in this section 
hold. Consider the initial value problem 

A = AX('{t,X), 

AiO) = h, ^^-^^ 

where X'^{t,X) defined by j^.^] ) is sufficiently CS and ^^-periodic in t for A > small. Then, 

there exists a sufficiently CS solution A{t,X) = e'^^^'^^B{t, X) of the initial value problem J^.dp , 
where X{X) G so(3) and B{t,X) is a -^^-periodic function such that 5(0, A) = I3. 



We call B{t, A) the periodic part of $(i, ua, A) associated to X{X). Using Lemma we get the 
following result: 

Theorem 4.4 (Hopf Bifurcation Theorem for Rotating Waves on a Sphere). Suppose that 

the assumptions made in this section hold. Then, there exists a sufficiently CS branch ^{t,ux, X) 

of solutions for the reaction- diffusion system J^. j| ) such that ^(t, uo,0) = c^'^^uq and for A > 

small, ^{t,ux, X) is either an orbitally stable modulated rotating wave with a primary frequency 
> 

vector X{X) and the secondary frequency uJx, or an orbitally stable periodic solution with the period 

A slightly modified version of this theorem was proved by |4H IHH] . They used the Liapunov- 
Schmidt reduction on scales of Banach spaces. We use the equivariant manifold reduction and 
Lemma 14.31 

The bifurcation diagram is illustrated in Figure ETTl 
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amplitude 



stable MRW 



stable RW 



unstable RW 



RW-rotating wave 

MRW -modulated rotating wave 



Figure 4.1: Bifurcation Diagram for Nonresonant Hopf Bifurcation 



5 The Primary Frequency Vector Formula for ^(t,ux,X) 

Let Xq = j^Xq. There exist Xi,X2 G so(3) such that the set {Xq, X i, X2} is an orthonormal 
basis in satisfying 

XlxXi = X2, XixX2 = Xl and X2 x = Xi. 



Since so(3) is isomorphic to M'^ and [X, y] = X x Y , {Xq, Xi, X2} is a basis of the Lie algebra 
so(3) such that 

[Xi, Xi] = X2, [Xi, X2] = Xi and [X2, X^] = Xi. 
Let us consider the reaction-diffusion system (j4.1|) . 

Throughout this section, we will suppose that uq G is a relative equilibrium that is not an 
equilibrium for (|4.1|) at A = and such that the stabilizer of uq is E„q = 13. Also, let $(t,no,0) = 



Suppose that Theorem 14.41 holds. Recall that a primary frequency vector X(A) of ^{t,ux, X) is 
given by ^(]^, A) = e for A > small. 

5.1 Existence of a sufficiently CS branch X{\) 

The theorems stated in this subsection can be proved without using BCH formula in so(3) |^. 
Let us define the smooth following function 



qm) 



Y 



if y , j^^Xq are in the same hemisphere or 1" = O3 



(— ^ + l)y ifiyyy, 7^X0 are in different hemispheres and Y ^ 



(5.1) 



\y\ 



where points and 27r are identified, and \Y\ is an equivalence class of D (see Appendix IXI for 
definition of D). Then e'^^l^D = el^l for any [Y\ G D. 
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Theorem 5.1. Suppose that the hypotheses of Theorem are satisfied. Then, there exists 
a unique sufficiently CS branch X[\) such that |^X(A) G (l{T)) for A > 0, X{\) is a primary 



frequency vector of ^{t,u\, X) for A > small and 



2n 

^bif I 



,no,0) 



X(0)i 



no- 



The branch ^(A) does not have the property that X(0) = Xq. The next corollary shows that 
we can construct a branch of primary frequency vectors X^ (A) such that X^ (0) = Xq. 

Corollary 5.2. |^ Suppose that the hypotheses of Theorem \4-4\ ^'"c satisfied. Then, there exists 

a branch X^ (A) for A > such that X-^(A) is a primary frequency vector of ^{t,ux, X) for A > 
small and X^ (0) = Xq. Moreover, 

1. if \Xq\ ^ kuobif for all A; G Z, then the branch X^ {X) for X> is sufficiently CS and 

(A) = \\Xq\ + 0{^/X)\ Xl + 0{^/X)Xl + 0{y/X)X2 for small A > 0. (5.2) 

(The branch ^{t,ux, A) contains only modulated rotating waves for A > small). 

2. if \Xq\ = kiOhif for some k Z, k 0, then |X'^(A)| for X >0 is continuous and 

Xf{X) = \Xo\ + 0{Xi) for small A > 0. (5.3) 

Therefore, we can not have frequency -locking phenomena for the modulated rotating waves 
^{t,u\, X). If |Xo| = kujhif for some /c G Z, 7^ 0, it may be possible that the branch X-^(A) be 
discontinuous at A = or/and at any A > such that |^(A)| = 0. 

5.2 Construction of the branches X(X) and X^(A) 

In the rest of the section we give a way of constructing X{X) and X-f (A) for A > using the BCH 
formula in so(3) presented in Theorem IA.2I Theorems 15.31 15.41 and 15.51 are used to construct the 
branch ^(A). Recall that a primary frequency vector X{X) of $(t, n^, A) is given by A{j^^,X) = 

for A > small, we show that there exists a sufficiently CS function Z{t, A) G q{D) such 
that A{t, A) = e^(*'^) for at least t G 



n ^ 



and A > small. 



Theorem 5.3. Suppose that the hypotheses of Theorem \4-4\ o.re satisfied. Let us consider the 
following initial value problem 



Z 
Z{0) 



1 



COS ^ 

2sinlfi|Z| 



^^(i,A), 



(5.4) 



Then, 



1. there exists a positive integer n independent of X such that the initial value problem \5.4\ has 



a unique sufficiently CS solution Zi{t, A) on t E [0, 



2lT 



and A > small; 
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2. 



Zi{t,\)= (^\Xo\t + VxJ^ xoi{s)ds + Xxo3{t,X)^ X^ + ^zi{t,\)Xi + ^Z2{t,\)X2 (5.5) 

for any t G [0, ^^J^ | ] and A > small; 

The integer n that appears in Theorem 15 . 31 can be found and, in general, it is not 1. Therefore, 
in general, the function Zi{t, A) is not defined on the entire interval [0, ^^], for A > small. 

Corollary 5.4. Suppose that the hypotheses of Theorem \4-4\ satisfied and the positive integer 
n is the one obtained in Theorem \5.^A Then, for any i = 1, 2, . . . , n — \ the initial value problem 



Z 

2-K 



^3+2^+(^-2sinM|z| 



^ ' (5.6) 



has a unique sufficiently CS solution Zi^i(t, A) on [i ^^J^ | , [i + 1) ^^J^| ] and A > small. 

Combining Theorem 15.31 and Corollary 15.41 we obtain the following result: 

Theorem 5.5. Suppose that the hypotheses of Theorem \4-4\ ^'"e satisfied. Then, there exists a 
sufficiently CS function Z{t,X) £ q{D) such that A{t,X) = e^^*''^^ for t £ [0,cxd) and A > small, 
where A{t, A) is the solution of the initial value problem i4.^ . In fact, the function Z(t, A) satisfies 
the following initial value problem for t £ [0, oo) and A > small : 



Z(0) = O3 



\Zr 2sin^|Z| 



XG{t,X)i modlTT), ^^^^ 



Remark 5.6. In fact, the function Z{t,X) obtained in Theorem \5.5\ is constructed as follows: 
1. On the interval [0, —I.' 

Z{t, A) = Z\t, A) = q{BCH{Zi{—- A), ^2(2^—, A), . . . , A))) (5.8) 

for all t £\{i — 1) ?^ I , i ?^ , 1 and A > small, where i = 1, 2, . . . , n and Zj is the solution 



of the initial value problem l\5.f)]) on the interval [{i — 1) 



27r • 27r 



2tt (■ I 1 \ 27r 



2. On the interval [^j^^^) (^ + /''^ integer i > 1 and A > small : 

?7r 27r 
Z{t, A) = Z\t, A) = q{BCH{Z\- — X),Z''\t - . — A))). (5.9) 

Using the function Z(t,X) obtained in Theorem 15.51 we construct the branches X(X) (A > 0) 
as follows: 
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Proposition 5.7 (Primary Frequency Vectors Associated to ^{t,ux, X)). Suppose that the 
hypotheses of Theorem are satisfied. Then, the branch X{X) for A > from Theorem \5.1\ is 
defined by: 

X(X) = A) for small A > 0, (5.10) 

zvr I^aI 

where Z{t,X) is given in Theorem \5.5\ 

Corollary 5.8 (Primary Frequency Vectors Associated to ^{t,u\, X)). Then, the branch 
Xf (A) for A > from Corollarv \5.^ is defined by: 

1. if \Xo\ / kujuf for all k e Z, then we define Xf{X) = '^^ff X(A) /or A > small. 

2. if \Xq\ = kujbif for some G Z, /c / 0, then we define: 
Xf{0)=Xo; 

if \X{X)\ / 0, then we define Xf{X) = ^-^^^^}^^f^ X {X) for X > small; 

if \X{X)\ = 0, then we define X^{X) = X{X) + k\ujx\ QW, where Q(A) G so{3), \Q{X)\ = 1 
for A > small. 

Using Lemma lOl for A > small, we can associate A{t, A) = e^('^)*5(t, A) = e^^^^'^^B^ (t, A) to 
each modulated rotating wave or periodic solution of period |^ obtained by Theorem 14.41 where 
B{t, A), B^' (t, A) are |^-periodic in t for A > small, as well as B{0, X) = B^ (0, A) = Is. Also, let 
B{t,0) and Bf{t,0) be such that A{t,0) = e^^^'^^B{t,0) = e^-^W*S/(t, 0), that is Bf{t,0) = I3. 

Proposition 5.9 (Periodic Parts Associated to ^{t,ux, X)). Suppose that the hypotheses of 
Theorem^^ are satisfied. Then, 

1. if \Xq\ 7^ kujhif for all k £ Z, we have 

(t. A) = e^'*"^*'^) for t G [0, 00) and X > small, (5.11) 
where Y{t, A) is -^^^-periodic in t and sufficiently CS for t G [0, 00) and A > small. 

2. if = kujhij for some k £ Z, k ^ 0, we have 

B{t, A) = e^o*+^-f^(*'^) for t G [0, 00) and X>0 small, (5.12) 

where H{t, A) is sufficiently CS for t G [0, 00) and A > small, and e 1^ = I3 

for A > small. 

Remark 5.10. In fact, Bf{t,X) = ei^^^^^i-xs {\%z{t,m (respectively B{t,X) = ei(BCHi-X{x)t,z{t,\))) j 
for t G [0,00) and A > small, where Z(t, A) is given in Theorem \5.fA 

Using Theorems 15 . 71 and we get the following result: 

Theorem 5.11. Suppose that the hypotheses of Theorem \4.4\ satisfied. Then, 
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1. if \Xq\ 7^ kiVbif for all A; € we have 

*(t,nA,A) = e[(l^ol+«(^))^o+^(^)^i+«(^)^^]*e^^(*'")^(VAr(t, A)), 
where Y{t, A) and r{t, A) are sufficiently CS and ^j^-periodic in t for A > small. 

2. if = /cwftij for some k £ Z, k ^ 0, we have 

^it,ux,X) = e[o(v^)^o+«(v^)^i+0(v^)^^]*e^»*+^^(*'^)M/(VAr(t,A)), 

where r{t, A) and H{t, A) are sufficiently CS and r{t, A) is ^j^-periodic in t for A > small, 

Xot^+VAH(t^,A) 
as well as e 1"^^! 1"^^! = Is for A > small. 

In fact, we have that the initial value problem ()4.2|) is equivalent with the following initial value 
problem 



Z 



^3 + ^^+(^-^)^^ 



XG(g,A)( mod 27r), 

(5.13) 



4 = XN{q,X), 

Z{0) = 'Os- 

We could have considered that the differential equation q = X]\f{q,X) has a |^-periodic solution 

q{t, A) = 0(\/A) instead of assuming that a supercritical Hopf bifurcation takes places at g = for 
A = 0. All the results remain valid. 



6 Examples 

Throughout this section, we let e = \/A. Here we present some examples of functions X^{t, A) for 
which we can get the closed form solution A(t,X) of the initial value problem 1)4. 5p . They obey 
Theorem 15.111 

Example 6.1 ( Hopf bifurcation to modulated rotating waves). Let g{t,X) be a sufficiently 
smooth periodic function of period -, — . such that g(0, A) = for A > small. For any t G [0, oo) 

and A > small, we define 

A) = (26Xi + 2eX2 + 2eX'o)g{t, A) + e-(2^^i+2eX,+eXl)<,{t,A)(^^ ^ ^^^)g(2.X,+2.X,+2.Xi),(t,A)_ 

It is clear that X'^(t,A) is a sufficiently CS, -, — i -periodic function such that X'~'(t,0) = Xq. 

|'^M/+A| 

Then, the initial value problem i4.^ has the solution 

A{t,X) = e(-^o+eXi)tg(2eXi+2eX2+2eXl)g(t,A)_ 

Clearly, Xq + eXi is not orthogonal to Xq. 

Example 6.2 (Example 1. Resonant drift phenomena for modulated rotating waves). 

Let uJbif = \Xo\ and g{t,X) be a sufficiently smooth periodic function of period — , such that 
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£,(0, A) = /or A > small. 

For any t £ [0, oo) and A > small, we define 



X^it,X) = {Xo + eX2)\^^^Ml^ + Xg{t,X)] +ee Jx 

{Xo+eX2) ( t+\g{t,X) ] 

.e V\/l^oP+A J, (6.1) 

is clear that A) is a sufficiently CS, -i — i - periodic function such that 0) = Xq. 



Then, the initial value problem ^'^^ ihe solution 



^(t, A) = e^^i*e VVl^oP+A ' (6.2) 

Clearly, eXi is orthogonal to Xq. 

Example 6.3 (Example 2. Resonant drift phenomena for modulated rotating waves). 

Let uji,if = \Xq\ and g(t, X) be a sufficiently smooth periodic function of period i — i such that 

g{0, X)=0 for X>0 small. 

For any t E [0, oo) and A > small, we define 



X°(t, A) = (.Y„ + tX^) I ' '•' + \g{t. A) I + ee 



-{Xo+eX2)\ '"'"■^t, t+\g{t,\) 



^ol' + A 



{Xo+eX2)\ ' f+Ag(f,A) 



is c/ear t/iat X'~'(t,X) is a sufficiently CS, ^ — i -periodic function such that X'~'(t,0) = Xq. 

|'^M/+A| 



■{Xo + Xi)e VVl^o|2+A y. (6.3) 

+a[ 

Then, the initial value problem i4-5\ ) has the solution 

^(i, A) = e^(^o+^i)*e VVl^ol^+A y. (6.4) 

Clearly, e(Xo + Xi) is noi orthogonal to Xq. 

Example 6.4 (Example 3. Resonant drift phenomena for modulated rotating waves). 

Let uJbif = ^ \Xo\ for /c G Z, A; 7^ and g{t, X, fi) be a sufficiently smooth periodic function of period 
1 — I such that g(0, A, /i) =0 for A > 0, small. 

For any t E [0, 00), A > small and small, we define 
X''{t,X,^,) = {Xo + ^iX^) | Ab^±^ + A5(t, A, /.) 



+ e'=e VVi^ol^+M^ ; ((e _ ^)Xo + Xi + X2)e 



^{Xo+^iXl)[ -!^]^^i^^t+Xg{t,\,^i)] (Xo+AtXi)( 4t^==t+A!?(t,A,M) 



(6.5) 
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It is clear that X'^(t, A, /i) is a sufficiently CSx, — , -periodic function such that X'~^(t, 0, /i) 

\^bif+-^\ 

Xq + fj,Xi . Then, the initial value problem i4-5\ ) for two parameters X, /i has the solution 



Clearly, for ji = \/A, we have orthogonality, that is e^{Xi + X2) is orthogonal to Xq. 

Example 6.5 (Nonuniform rigidly rotation). Let ajhij = \Xq\ and g{t,X) be a sufficiently 
smooth periodic function of period -, — . such that g(0, A) = for A > small. 

For any t € [0, 00) and A > small, we define 

X^{t, A) = (Xo + eXi)(l + eg{t, A)). (6.7) 

It is clear that X'^{t,\) is a sufficiently CS, — , -periodic function such that X'^(t,0) = Xq. 
Then, the initial value problem i4.5\ ) has the solution 

^(t,A) = e(^"+^^i)(*+^?(*'^)). (6.8) 



This clearly represents a non-uniformly rigid rotation about the line containing the vector Xq + eXi . 
This is a degenerate situation for the case of one parameter. 

Also, the case of a |^-periodic solution is a degenerate situation for one parameter. Therefore, 
we do not present it here. 

These examples will be used in Section |S| As we have seen in Examples 16.21 and 16.31 if we have a 
single parameter A and \Xq\ = kuj^if for some k £ Z, we can have three different types of solutions 

for the initial value problem (|4.5|) . In Example (|6. 4 j) . we have two parameters A > and /u, as well 

> 

as |Xo| = kiObif for some A; G Z and we have a branch fx = /u(A) for which X{X, fi{X)) is orthogonal 
to Xq. In next section, we prove that this is generically valid. 

7 Resonant Drift Phenomena for Modulated Rotating Waves on 



We consider the following reaction-diffusion system 

rill, „ 

-{t, x) = DAsu{t, x) + F{u{t, x), A, /x) on rS^, (7.1) 



dt 

( di ... \ 



with dj > for 



where n = (ui, ^2, . . . , uat) : M x rS^ with N > I, D = 

V ... dN J 

i = 1,2, iV are the diffusion coefficients and F = (Fi, F2, . . . , Fn) : x M x M ^ are 
sufficiently smooth functions such that F(0, A,/i) = for |A|, small. 

We study the reaction-diffusion system 1)7. 1|) on the function space Y defined in Section |21 Let 
$(t, u, A, fi) be the 50(3)-equivariant sufficiently smooth local semiflow defined as in Section|2l Let 
uq £ Y° be a relative equilibrium that is not an equilibrium for (|7.1|) at (A, ^) = (0, 0) and such 
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that the stabihzer of uq is = I3. consider L the Unearization of the right-hand side of (|7.1() 
with respect to the rotating wave no,0, 0) = c^^^uq at (A,/i) = (0,0) in the co-rotating frame, 
that is 

L = DAs + DuF{uo,0,0)-Xo. 

Suppose that: 

1. a{L)ri{z £ C I Re (z) > 0} is a spectral set with spectral projection P*, and dim(i?(P*)) < 00; 

2. the semigroup e^* satisfies |e^*|_R{i_p^) | < Ce~^°* for some f3o > and C > 0. 
Theorems IB . 1 1 with parameters A, /i can be applied. 

There exist sufficiently smooth functions Xq '■ x M x R ^ so(3) and X^- : 14 x M x M ^ such 
that any solution of 

A = AXG{q,X,^i), .^2) 
q = XAr(g,A,^), 

on SO{3) X Ki, corresponds to a solution of the reaction-diffusion system H4.1() on M^^{X,n) under 
the diffeomorphic identification for |A| and |/u| small. Also, Xg{0,0,0) = Xq, XAr(0,0, 0) = 
and (j(D„Xjv(0, 0, 0)) = o"(Q*i|v;)i where Q^: is the projection onto K along Tuq{SO{3)uo)- Let 
dimi?(P*) = 5. Suppose that a supercritical Hopf bifurcation with eigenvalues ziziuJuf takes place 
in the second differential equation of 1)7. 2(1 in 14 at g = for (A, /u) = (0,0), that is: 

1. X7v(0,0,0) = 0; 

2. DqXiy{0, 0, 0) has eigenvalues iiujbif, without loss of generality, we assume that Xn{0, A, /i) = 
for |A|, \iJ,\ small; 

3. DgXiy{0, X, fi) has the eigenvalues a(A, ^u) ±i{u!hif +/3(A, ^) with a;(0, 0) = /3(0,0) = 0) and 
aA(0, 0) > 0. This implies that a^XH^lj)-, fj) = for some sufficiently smooth curve A = XnifJ-) 
with Xh{0) = 0. This curve represents the Hopf points and without loss of generality, we 
suppose that A/f(/i) = for small; 

4. the branch of periodic solutions q{t, A, fi) bifurcating from q = generically satisfies q{t, A, /x) = 
y(A)r(t,A,/x). 

For A > small and Ifxl small, let T(A, fi) = , , be the period of the periodic solution q(t, A, fi) 

near q = 0, that appears due to the supercritical Hopf bifurcation, where wa,^ = "^m/ + 0(m) + 
As(A,/i) for A > small and small. 

By Theorem 14. 4( there exists a sufficiently CS branch #(t, liA,/^, A, fi) for A > small and |^| small 
such that for |//| small, ^{t,uo,fi,0, fi) = c-^^^^'^'^^^uq^^ and such that for A > small and |^| small, 
^{tjUx X, ^) is an orbitally stable modulated rotating wave or periodic solution of period , 1 . 

Let X(A, ^) be the sufficiently (75 branch such that for small, e '"''•''I =e ''^''•''l , and 

> 

for A > small and |//| small, X{X, fi) is a primary frequency vector corresponding to ^{t, ma,/^; A, /i). 
We have the following two results: 

Theorem 7.1 (Resonance Case). Suppose the previous assumptions hold and that \Xq\ = kio^if 
for some k £ Z, k ^ 0. Let 

^g(9,A,^) = Xo{q,X,fi)XQ + xi{q,X,fi)Xi + X2{q, X, fi)X2. 
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X 



stable MRW 



X (n)=0 



stable RW 



represents MRW that are orthogonal to 
or periodic solutions 



Figure 7.1: Parameter Space for Resonant Hopf Bifurcation 

Then, if (xq) ^{0,0,0) / k{uJo^fj_)' \fj_=o, there exists a sufficiently CS curve ^ = ^(A) for X > small 

such that /i(0) = 0, and there exists a sufficiently CS branch of orbitally stable modulated rotating 

> 

waves ma,/^(a)) /^(-^)) for A > small, with a primary frequency vector X(\, //(A)) orthogonal 
to Jto. 

Corollary 7.2. Suppose that Theorem \7.1\ hold. Then, 

k 

1. the branch X{X,iJ,(\)) is sufficiently CS for A > small and X(A, /u(A)) = 0(A)2)Xi + 
0(At)X2 for X>0 small. 

2. fi(t,A,/i(A)) = e(0{Vx)x^o+o{^)x,+o{^/\)X2)t^Xot+VxHit,x)^ ^ ^ [0,cx)) and X > small, 
where H{t, A) is a sufficiently CS function such that e = for A > small. 

3. there exists a branch (A,/i(A)) with A > that is discontinuous at X = 0, but \X-f (A, /i(A))| 
is continuous for A > small and |X-'^(A, ^(A)) | = |Xo| + 0{Xi) for A > small. 

The branch ^ = /i(A) for A > smaU can be found using the BCH formula in so(3). The 
parameter space in the case of the resonant Hopf bifurcation is illustrated in Figure H7.1() . 

The results in Sections |1J [3 and [7| show that a rotating wave on a sphere generically under- 
goes a transition to quasi-periodic meandering (modulated rotating wave) at a Hopf bifurcation. 
Furthermore, the primary frequency vector of the modulated rotating wave is determined by both 
the critical Hopf eigenvalues and the frequency of the rotating wave undergoing the bifurcation. In 
particular, resonances between the critical Hopf eigenvalue and the frequency of the rotating wave 
undergoing the bifurcation lead to orthogonal drift. While we have concentrated on the case of 
Hopf bifurcation, it should be mentioned that we would obtain similar results for periodic forcing 
of a rotating wave: equation (|4.5)) represents the dynamical equations of a periodically forced ro- 
tating wave. In a spherical heart, for example, one could alter the dynamics of a rotating spiral 
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wave simply by applying a (weak) periodic forcing (e.g. pacing), thereby inducing meandering. 
By an appropriate choice of the forcing frequency, one could then control the direction of the 
primary frequency vector (i.e. steer the spiral tip to another location). Of course, this assumes 
perfect spherical symmetry. The case when geometrical imperfections appear can also influence the 
dynamics of rotating waves and it will be presented somewhere else. 



8 Some Numerical Results 

Theorem 15.111 shows that, generically, by a supercritical Hopf bifurcation of a rotating wave, not 
only we get modulated rotating waves (that is a quasi-periodic tip motion) , but also we get a quasi- 
periodic meandering tip motion. It is possible that this meandering motion is not of epicycle-type. 
We illustrate this using the examples presented in Section 1^1 

Let Lt = 



Rz{0) 




We have seen that the study of the reaction-diffusion system 1)4. 1|) on Y" reduces to the study 
of the finite-dimensional system (|4.2j) on the center manifold M^"(A) of the relative equilibrium 
50(3)^0, with *(t. Mo, 0) = e^o*tio. Let us denote ^^(g, A) = F^(g, A)L^. + F?'(g, A)Ly + F^(g, A)L^ 
for any g G K and |A| small. 

If we parameterize 50(3) by Euler angles, that is ^ = Rz{'4')Rx{d)Rz{4')-, where (f) G [0,27r), 
G [0,7r], ^ G [0, 27r), the finite dimensional system (|4.2)) becomes 



= F^(g,A) -cot0[FS^(g,A)cos</) + F^(g,A)sin,/.] , 
e = -F^(g,A)sin0 + F-(g,A)cos0, 
^ = ii^[i^n9,A)cos</. + F-(g,A)sin(/>], 
q = XN{q,X)- 

We consider the following initial value problem associated with the system (|8.1j) : 

= F^{q,X) -cot 9[Fy{q,X) cos ^ + F''{q,X) sin 

9 = -Fy {q,X) sin (P + F''{q,X) cos (p, 

^ = ^[Fyiq,X)cos<P + F-{q,X)sin<p], 

m = 0, (8.2) 

9{0) / 0, 

V'(O) = 0, 

q = XN{q,X). 

We choose 6'(0) / near in (lO) . 

We consider K ~ C and a supercritical Hopf bifurcation takes place in q = Xiy{q,X) at q = 
for A = with eigenvalues ztzw^jj-. Let q{t,X) be the periodic solution of the second differential 
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equation in (|B.1|) . that appears by a supercritical Hopf bifurcation for A > small. Let T{\) = 
be its period, where ujx = ujuj + 0(A) for A > small. Let q{t,Q) = for any t G R. If we 
substitute g(t, A) in Xq^q, A) and in the first three equations of the system (|8.2() . we get X^(t,A) = 
Xciqit, A), A) = F-{q{t, A), A)L,. + Fy{q{t, A), X)Ly + F-{q{t, A), A)L, = |Xo| + ^/\H{t, A) for 
A > smah and t £ R. If we write F^^(t,A) = F^(g(t, A), A), Fyy{t,X) = Fy{q{t,X),X) and 
F^^(t, A) = F^(g(t, A), A), we get that X^{t, A) = F^^(t, X)L^ + Fyy{t, X)Ly + F''{t, X)L^ is 
periodic. 

The following initial value problem is obtained 

= F^^(t,A) -cote[F2'f(t,A)cos0 + F^^(t,A)sin(/.], 
= -Fyy {t, X) sin (p + F'''={t,X) cos (f), 

^ = ii^[i^^"(t,A)cos0 + F--(t,A)sin(/>], , . 

m =0, ^ • ^ 

9{0) / 0, 
-(/-(O) = 0. 

We choose 6'(0) / near in (|01) . 

We use a Maple program that integrates numerically the system ()8.3() and finds numerically a 
primary frequency vector (|llj, if its norm is not or vr) of the modulated rotating wave that 
appears by a supercritical Hopf bifurcation, as discussed in Section \^ Then, for a choice of the 
point xq S rS^ near j^Xo, we represent the tip motion 

xupi^it, ux, A)) = ^(0, X)-^A{t, X)xo on rS^ (8.4) 

for A > small. Recall that ^{t,ux,X) = ^(0, A)"M(t, A)*(g(t, A)). 
We consider the following cases: 

Case 1. X{t,X) is given in Example ()6.1|) from Sectional Xq = Lz,Xi = Lx,X2 = Ly] totif = 20, 

|Xo| = 2, g{t,X) = sm{{uJbif + X)t), r = 3, 9(0) = 0.01, xq = 0.92 | . We get Figures lO 

V 2.85 

toOfor A = 0.01, 0.05. 
Since l-'^^ol / kuJbif for any k £ Z, we have the nonresonant case. On Figures ISH and l8?2l it is 
visualized the tip motion given by 1)8. 4|) of the modulated rotating waves $(t, ux, A) obtained 
by a supercritical Hopf bifurcation that takes place in q = XN{q, X) at q = for A = 0. On 
Figure IHTTI we consider A = 0.01. On Figure IH?^ we consider A = 0.05. 

On the second of Figure WJ\ we plot the points xtip{^{i^,ux, X)) for some i = 0, 1, 2,. . . 
and we can see that they are points of a circle on the sphere rS^ with the center on the line 
having the direction of the primary frequency vector of $(t, ux, A), where A = 0.01. 
The grey line is the line containing the frequency of the rotating wave undergoing the Hopf 
bifurcation, Xq. The black lines are the lines corresponding to the primary frequency vector 
associated to the modulated rotating wave. One of them is computed numerically and the 
other is the exact one. We can see that they are very close. 

Case 2. X{t,X) is given in Example (|6.2|) from Sectional X^ = Lz,Xi = Lx,X2 = Ly] uj^if = \Xq\ = 



20, g{t, X) = sin[{ojbif + X)t), r = 3, 61(0) = 0.02, xq = \ 0.92 

2.85 
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We get Figures lO to lO for A = 0.05, 0.1. 

Since |Xo| = kujbif for /c = 1 G Z, we have the resonant case. On Figures EHl and E31 it is 
visuahzed the tip motion given by 1)8. 4|) of the modulated rotating waves ux, A) obtained 
by a supercritical Hopf bifurcation that takes place in q = X]y{q, A) at g = for A = 0. On 
Figure IH^ we consider A = 0.05. On Figure IHUl we consider A = 0.1. We can see that the 
primary frequency vectors associated to the modulated rotating waves are orthogonal to Xq. 
We call this phenomenon resonant drift. 

On the second of Figure ESI we plot the points xup{^{i^,u\, X)) for some i = 0, 1, 2,... 
and we can see that they are points of a circle on the sphere rS^ with the center on the line 
having the direction of the primary frequency vector of u\, A), where A = 0.05. 
The grey line is the line containing the frequency of the rotating wave undergoing the Hopf 
bifurcation, Xq. The black line is the line corresponding to the primary frequency vector 
associated to the modulated rotating wave. 

Case 3. X{t,X) is given in Example (|6.3|) from SectionlHl X^ = Lz,Xi = Lx,X2 = Ly] Uf^f = \Xo\ = 

/ 0.44 

20, g{t, A) = sin((u;;,i/ + A)t), r = 3, 9{0) = 0.5, xq = 0.14 

V 2.96 

We get Figures EH to Eli for A = 0.05, 0.25. 
Since |Xo| = kuJbif for /c = 1 G Z, we have the resonant case. On Figures and IHl^ it is 
visualized the tip motion given by 1)8. 4|) of the modulated rotating waves ux, A) obtained 
by a supercritical Hopf bifurcation that takes place in q = Xjv(g, A) at g = for A = 0. On 
Figure we consider A = 0.05. On Figure EUl we consider A = 0.25. We can see that the 
primary frequency vectors associated to the modulated rotating waves are not orthogonal to 
Xq. This can happen in the resonant case if we consider only one parameter A. 
On the second of Figure we plot the points xtip{^{i^, ux, X))ior some i = 0, 1, 2,... 
and we can see that they are points of a circle on the sphere rS^ with the center on the line 
having the direction of the primary frequency vector of ux, A), where A = 0.05. 
The grey line is the line containing the frequency of the rotating wave undergoing the Hopf 
bifurcation, Xq. The black line is the line corresponding to the primary frequency vector 
associated to the modulated rotating wave. 
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20T)'"o.oi,0.01)),i — 1, ... ,5 

Figure 8.1: Case 1, A = 0.01 

The nonuniformly rigid rotation shown in Example 16.41 is not generically. Therefore, we do not 
present it here. 

9 Proofs of Theorems 

Proof of Lemma \4-^ Let us omit the parameter A. Since SO{3) is a compact manifold and X*^ 
is a sufficiently CS function, the initial value problem ()4.6|) has a unique sufficiently CS solution 
that is globally defined. Let A*{t) be this solution. Since the exponential map exp: so(3) — > S0{3) 
is surjective, there exists a matrix X G so(3) such that A*{T) = e^'^ . 
Let us make the following change of variable B = e~^^A. 
Then the first equation of 1)4. 6|1 becomes: 

B = -e-^^XA{t) + e-^*A{t), (9.1) 

or since A = AX^{t), we get 

B = -e-^^XA{t) + e-^^AX^{t), (9.2) 
or since A = e^^B it follows that 

B = -XB{t) + B{t)X^{t). (9.3) 

The initial condition ^4(0) = Is becomes 5(0) = /s- 
Let us consider now the initial value problem 

B = -XB + BX^it), 

B{0) = h. ^^-^^ 
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Xtip(*(t,Mo.05,0.05)) 



Figure 8.2: Case 1, A = 0.05 

The initial value problem H9.4|) has a unique global solution since S0{3) is a compact manifold. 
Since A*{t) is the solution of the initial value problem (|4.()j) . it follows that B*(t) = e~-^*A*{t) is a 
solution of (EiJ and B*{T) = I3. 
Let us define 

C*{t) = B*{t + T). 

We will check that C* is another solution of (|9.4|1 . We have C*(0) = i?*(r) = and for any 
t G [0,00), 

dt = + = + = + T) + B*{t + T)X^{t + T) ^g^^^ 

= -XB*{t + T)+ B*{t + T)X^{t) = -XC*{t) + C*it)X'^{t), 

where we have used the fact that X^ is a T-periodic function. 
Therefore, B*{t) = C*{t) for t G [0,oo), that is, B* is T-periodic. 
So, = e^^B*{t) for t G [0,oo). 

Taking into account the parameter A > 0, since X^{t, A) is sufficiently C5 for t G [0, 00) and A > 0, 
we have that A{t, A) is sufficiently CS for t G [0, 00) and A > small. If we work with A > 0, then 
A{t, A) is sufficiently C5 for t G [0, co) and A > small. □ 

Proof of Theorem \4-4\ We apply Theorems IB. II Then, it is sufficient to study the differential 
equations (|4.2j) . 

The supercritical Hopf bifurcation in V* for the second differential equation in (|4.2|) at g = for 
A = implies that there exists a unique sufficiently CS branch of periodic solutions q{t, X) for A > 
near q = 0. 

Recall that we have denoted by T(A) = the period of q{t,X), where u>x = ^Uf + 0(A) is a 
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xtipi^it, uo.05, 0.05)) xtipi^ii^^,uo.o5,0-05)),i = 1, . . . , 5 

Figure 8.3: Case 2, A = 0.05 
sufficiently smooth function for A > small. 

If we substitute q{t, A) in the ffi'st of the differential equations ()B.1|1 . we get the differential equation 

A = AX^{t,X), (9.6) 

where we recall that the function X^{t,X) is sufficiently CS and |^-periodic for t G [0,cxd) and 
A > smah (and it is defined by 

We consider the following initial value problem on SO {3): 

^(0) = Is, ^'-^^ 

because any solution of the differential equation (|9.6() with ^(0) = Aq, where S 50(3) is given 
by AQA*{t, A) and A*{t, A) is the solution of the initial value problem (|9.7j) . 
For A > small, by applying Lemma 14.31 we get 

A{t, A) = e^(^)*S(t, A) for t e [0, oo) and A > smah , 

where B{t, A) has period T(A) and -8(0, A) = I3. 

Consequently, by using Theorems IB. II there exists a solution ^{t,ux, X) of the reaction-diffusion 
system (|4.1() of the form: 

*(t, ux, A) = A{t, A)^'(g(t, A)) = e^(^)*S(t, A)^'(g(t, A)), where t G [0, 00) and A > smah , 

with ^'(g(0,A)) = ux and ^' appears in the statement of Theorem IB. II fsee Appendix IB)) . If we 
write 

q,{t,X)=B{t,X)^{q{t,X)), 



25 



Xtip(*(t,Mo.l,0.1)) 

Figure 8.4: Case 2, A = 0.1 

then the function qi{t, A) is T(A)-periodic in t, since B{t, A) and q{t, A) are T(A)-periodic in t. Also, 
we get 

*(t, ux, A) = e^W*^i(t, A) and *(r(A), ua, A) = e^W^^tXA- 

We check that, at least generically, ^{t,ux,X} ^ S0{3)ux for any t G (0,T(A)) and A > small. 

> 

We prove that X{X) is a primary frequency vector of n^, A). 

Suppose that there exists a Ti(A) G (0,T(A)) such that $(Ti(A), ua. A) G S0{3)ux and for t G 
(0,Ti(A)) we have ^{t,ux,X) ^ S0{3)ux. Let r(A) = lTi{X)+r, where r G [0,ri(A)) audi G Z. Let 
*(ri(A),iiA,A) = e^(-^)'^i(^)tiA, where y(A) G so(3). Using the definition and S0{3) -equivariance 
of the semiflow *, we get *(r(A),MA,A) = *(r, *(/ri(A), ua, A)) = e'^W^i W*(r, ua. A). 
Therefore, ^{r,ux,X) = e-'^(^)^i(^)e^(-^)'^(^)nA G S0{3)ux and r G [0,T(A)). It follows that r = 
and so T(A) = m(A). Thus, *(r(A),nA, A) = e^W^^UA = e^W^W^^A and e^W^W = e^W^W 
because = I3. 

Also, we get ^{t,ux,X) = e^W*^i(t, A) = e^W*e-^W*e^W*B(t, A)*(g(t, A)). 

Let Bi{t,X) = e"'^W*e^(^)*B(t,A). Si(t,A) is T(A)-periodic, because B{t,X) is T(A)-periodic and 

gX(A)T(A) ^ gr(A)T(A)^ rp^^g^ 

*(t,tXA,A) = e^W*i?i(t,A)^(g(t,A)). 

Then *(Ti(A),'UA,A) = e^W^^W^A = e^W^i WSi(Ti(A), A)^'(g(ri(A), A)), and it follows that 
^{q{Ti{X),X)) = {B{Ti{X),X))-W. 

We have that ux is close to uq and 50(3)7X0 H = {uo}, therefore, generically, S'0(3)nA n 

^(V*) = {ux}. Since *(9(Ti(A), A)) = {B{Ti{X), X))-^ux G "^{V^) n SO{3)ux, it follows that 
B{Ti{X),X)ux = ux and, since S,, = h, Bi(ri(A),A) = and *(g(ri(A), A)) = ux = ^(^7(0, A)). 
Because is a local diffeomorphism (see 0113) and q{t, A) = \/A?"(t, A), it results that q{Ti{X), A) = 
q{0,X). We know that q{t,X) is a solution of the differential equation q = Xiy{q,X). Therefore, 
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Figure 8.5: Case 3, A = 0.05 



it follows that q{t,X) is Ti(A)-periodic. But it has the period T(A). That is, T(A) = Ti(A) 
contradicting the fact that ri(A) G (0,r(A)). 

If |-^(A)| r(A) = 2k7r for some /c G Z, then ^{t,ux, A) is a periodic solution with the period T{X). 
Otherwise, ua, A) is a modulated rotating wave with a primary frequency vector X{\). 
By Proposition 14.21 the orbital stability of the modulated rotating wave or of the periodic solution 
obtained above is the same as the stability of the periodic solution q{t, A). 

Therefore, since the periodic solution q{t, A) for A > is stable, then the corresponding modulated 
rotating wave or periodic solution is orbitally stable as well. 

Since A{t, A) and q{t, A) are sufficiently CS and ^' is a local diffeomorphism, we get that ux, A) 
is sufficiently CS. The same result can be inferred if we take into account that ux = ^'(g(0. A)). □ 

Proof of Theorem \5.8[ We will use the following result (see |U| ) : 

Lemma 9.1. Let g{t, u) be a continuous function on an open connected set [oi, 6i) x M"*" C O C 
and such that the initial value problem for the scalar equation u = g{t, u) has a unique solution 
u{t) > on t G [ai,6i). ///: [ai,6i) x R" ^ M" is continuous and \\f{t,x)\\ < g{t,\\x\\) for 
t £ [ai,6i) and x £ M", then the solutions of x = f(t,x), ||x(ai)|| < n(ai) exists on [ai,6i) and 
Mt)\\ < u{t) forte [ai,6i). 

There exists a constant M > independent of A > small such that |X^(t,A)| < M for 
t G [0, oo) and A > small. For A > small and any \Z\ < vr, we have 

\z\ \ 



II/3II + 2 11^1 



1 



cos ■ 



\Z\ 



2sin^ 



Z\ 



X''{t,X) 



< 



+ 



\z\ 



+ 



COS ■ 



2sinlf |Z 




(9.8) 
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Xtip(*(t,Mo.25,0.25)) 



Figure 8.6: Case 3, A = 0.25 



h + -z + 



cos 



\z\ 



\Z\ 



2sin^ |Z| 



M 



V3 + ^V2|Z|+ 1^7 + 



cos 



\z\ 



\Z\ 



2sml|i |Z| 



2\Z\ 



< 



or, if we use the fact that the function x 



cos 



is decreasing on [0, 5 



2J' 



\Z\ 



\Z\^ 2sin^|Z| 



Z 



M 



V3 + ^V2|Z| + (2 + 2-1) 



< M{<o + \Z\ 



< 



The initial value problem 



a 

a(0) 



M(6 + o), 




(9.9) 



(9.10) 



(9.11) 



has solution on an maximal interval [0, bmax], where bmax is defined such that a{t) < vr on [0, bmax]', 
bmax is independent of A. 

Therefore, using Lemma l9. II with ai = 0, it follows that the sufficiently CS solution Zi{t, A) of the 
initial value problem (|5.4j) is defined for any t G [0, 6m,ax] and A > small, and that |Z(t, A)| < vr 
for any t G [0, bmax\ and A > small. 

If we choose a positive integer n > such that n > then the initial value problem (15.41) has 

Omax ' 

a solution Zi(t,X) defined and sufficiently CS on t G [0' ~?r^] and A > small. Since r(A) = 
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T(0) + 0(A), then we can choose n independent of A for A > smah, where T(0) = , i . 

This ends the proof of (1) of Theorem 15.31 

To prove the formula (|5.5j) . we write the Taylor expansion of 

Z*{t,€) = zoo{t,e)Xl + zii{t,e)Xi + Z22{t,e)X2, (9.12) 
zoo{t,e) = Zo{t) + Zi{t)e + e^Z2{t,e). (9.13) 

If we substitute (|9.12j) , (|9.13|) and (|4.4)) in the first of the differential equations given by (|5.4j) with 
A = e^, identify the orders of e in both sides and then take e = \/A, we get the formula 1)5. 5|) for 
Zi{t, A). For the sake of simplicity, we omit e. 
The result is 



+ 



cos 



\/^00+^ll+^22 



Zoo + ^11 + ^22 2 v^Zq^o + + z'i^ sin V4.+^^!i ±fI 



•(ZOO-'^O + Zll^l + Z22^2)^]2;o-'^0 + (^^iXi + 6X2X2. 

(9.14) 

Since Z* is sufficiently smooth, we get Z*{t, e) = Z*(t, 0) + eH{t, e) = |Xo| tX^ + eH{t, e). 
Therefore, zii(t, e) = ezi(t, e) and Z22(t, e) = ez2{t, e). In the right-hand side of (|9.14jl . the coefficient 
of X^ is xo{t, e) + e'^k{t, e) = |Xo| -|- exoi{t) + e^xo2(t, e) + e'^k{t, e) and thus zqo = \Xq\ + exoi(t) + 
e'^ki{t,e), where xq, xqi, xo2 are defined in (|4.4j) . Therefore, 

zoo(t, A) = |Xo| t + e xoi(s) ds + e^xo3(i, e). 
Jo 

This gives the formula (|5.5j) after we relabel zi{t, y/X), Z2{t, \/A) and xo3(i, "V^A) to zi(t. A), zi(t. A) 
and xo3(t, A). This ends the proof of (2) of Theorem 15.31 □ 

Proof of Corollary \5.4\ We have the same notations as in the proof of Theorem 15.31 From the 
proof of Theorem 15.31 we know that the initial value problem 

a = M(6 + a), 

a(0) = ^^-^^^ 

has the solution a defined on [0, bmax] with a{t) < vr on [0, bmax]- Let Iq = G [0, bmax]- Consider 
the initial value problem 

The solution to this problem satisfies a{t) < vr for any t £ [to, 2to]. 

We repeat the previous argument for to = i--^ for z = 2, 3, . . . , n — 1. From this and Lemma l9. 11 
with ai = to = i^^^ for i = 2, 3, . . . , n — 1, we get the conclusion of Corollary 15.41 □ 
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Proof of Theorem\FIB Recall that r(A) = |^ for A > small 



Part 1 On the interval [0,r(A)]: 

Consider the initial value problem 



A 

m 



AX«(t,A), 



(9.17) 



Let us make the change of variable A = near in the initial value problem (|9.17jl . 
Then, using Proposition IA.4I (2) we get the following initial value problem in Z: 



Z{0) 

En=o ^i-dzrz 

Z{0) 

En=o}^i^dzy^ 

Z{0) 



e^X^(t,A), 
XG{t,X), 

03, 



or, using again Proposition IA.3I (1) we get 

v^oo (-1)" yn 7 
Z^n=0 (n+1)!^ ^ 



Z(0) 



X«(t,A), 
O3, 



03. 



Using Proposition IA.3I (3). we have that 



n=0 



(n + 1)! 



/3 + cos|Z| - 1 ^ ^ \Z\ -sin|Z| ^2 



(|^^ will proved later. 

If we put into (nOTl) . we get 



Z(0) = 03- 



We will prove later that for any \Z\ < 2tt, 



\Zf \Z\^ J 2 



cos 



\z\ 



\Z\' 2sin^|Z|, 
Since we are looking for |Z| < vr < 27r, the system (|9.23p becomes 



Z' 



z 

Z(0) 



3^2 

03, 



\Zr 2sin^|Z| 



z' 



X«(t,A), 



(9.18) 
(9.19) 
(9.20) 



(9.21) 



(9.22) 



(9.23) 



(9.24) 



(9.25) 
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where we are looking for a solution Z such that |^| < vr on some interval [0,to(A)]. 
We now prove (|9.24j) . It is enough to check that 



/3 + 



cos|Z|-l |Z|-sin|Z| 2 



-Z + 



\z\- 



Z' 



COS 



\z\ 



\Z\' 2sin^|Z| 



Z' 



h- (9.26) 



Using Proposition Ol (3), that is Z^ = -\Zy Z and Z'^ = -\Zy Z"^, we have 



/3 + cos|Z| - 1 ^ ^ \Z\ -sin|Z| ^2 



\ZY 



cos 



\z\ 



\Z\' 2sin^|Z|, 



cos ^ \ 2 COS \ Z\ — 1 



iZf 2sinM|z| 



Z^ + 



cos \Z\ — 1 cos \Z\ — 1 I 1 



cos 



2|zr 

|Z| • 



|Z|-sin|Zl 2 , \Z\-sm\Z\ 



+ 



\ZY 

\Z\ -sin|Z| / 1 



l^r 2sin^|Z| 



^3 ^ I ^ I - Pill v2 _^ I ^ I - QUI I ^ I 



\Z\ 



2\Z\- 



COS 



I2I 



\Z\' 2sinM|z| 



Z^ 



h + 
1 



1 , cos|Z|-l cos|Z|-l , cos|Z| - 1 |Z|cos^ 1 ^ sin|Z| 



+ 



cos 



+ 



2 sin 



1^1 



2 ^ 2|Z| 



Z 



\z\ 



\Zf 2sin|^|Z| 



+ 



cos|Z|-l |Z|-sin|Z| |Z|-sin|Z| 1 



1^1 



|Z|-sin|Z| cos-^ 



1^1 2sinf|Z| 



2|Zr 



Z' 



+ 



1^1 izf 



■+ 



h + 



•sin^Mcos^ ^ sin|Z| 



\z\ 



sm 



1^1 



2|Z| 



Z + 



,2M 



fcosf 



M 

2 



2|zr 



sm 



1^1 



Z' 



- h- 
(9.27) 



The proof of (|9.22j) follows by using Proposition I A. 31 (3) : 



n=0 



/3 + 



-(2A; + 1)! 



k=0 



{2k + 2)\' 



6! 



'3! 5! ' ' 7! ' ' 



\z\' + ...]z 



Z' 



where we have used the Taylor expansions for sine and cosine. 



We apply Theorem lO to get ^(t,A) = e^i^*'^), for any f G [0,-^] and A > smah 



(9.28) 
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where Zi{t, A) is the sufficiently CS solution of (|5.4() . We make the change of variable B = 
(74(^^, A))~^^. We get that the solution of the initial value problem H9.17() is given by 

A{t,X) = e^'^^'^^B{t,X), 
where B{t, A) is the solution of the initial value problem 

B = BX^it,X), 

Using the above argument and Corollary 15.41 we have that B{t, A) = e^^^*''*') for any t € 
["rT^' ^"rT^l ^'^'^ A > small, where ^2(^1 A) is the sufficiently CS solution of ()5.6|) for i = 2. 
If we continue this, we get for A > small and any t £ [{i — j 

where Z^{t, A) is given by 

for any t G [{i — 1)-^^, i-^^^] and small A > for i = 1, 2, . . . , n, where Zi is the solution 

of the initial value problem (|5.6() on the interval [{i — 
Then A)] is defined by formula (|5.8|) . 
Also, from (|5.5|) in Theorem 15.31 we get 

Zi{t,X) = MXo|(t-(i-l)^) + VA y* ^^^^ xoi{s)ds + Xxi^{t,X)] 

+ VXz{ {t, X)Xi + ^/a4 (t, A)X2 (9.31) 
for any t £ [{i — i^^] and small A > for i = 1, 2, . . . , n. 

Because the BCH formula in so(3) is smooth from so(3) x so(3) into D (see Theorem IA.2p 
and Zi{t^ A) for i = 1, 2, ... , n are sufficiently CS, then [.Z'(t, A)] is sufficiently CS on [0, T(A)] 
and A > small. 

Part 2 On the interval [zT(A), {i + l)r(A)] for any integer i > 0. 
Since 

A{t + r(A), A) = e^W^W A(t, A) = A(T(A), A)A(t, A) = e^"(^W'^)e^"(*'^), 

it is easy to see that we can define for any t G [^(A), 2T(A)] and A > small, 

[Z\t, A)] = BCH{Z\T{X), A), Z\t - T(A), A)). 

Because the BCH formula in so(3) is smooth from so(3) x so(3) into D, Z^{t, A) is sufficiently 
CS and r(A) is sufficiently smooth, then [Z{t, A)] is sufficiently CS for any i G [r(A), 2r(A)] 
and for A > small. 
We then repeat the above argument. 
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It is clear that the function A)] is sufficiently CS, and A{t^\) = e^^^'^^ for t G [0,cxd) and 
A > small. Let us define Zi{t,X) = q{[Z{t, X)]). Therefore, Zi{t,X) is sufficiently CS and 

Thus, Zi{t, A) is the solution of the initial value problem (|5.7j) . 

□ 

Proof of Proposition \5. 7\ If we define 

X{X) = 7^Z{T{X),X) for A > small , (9.32) 

then ^(r(A),A) =e^W^W. 

Since Z{t, A) is sufficiently CS for t £ [0, oo) and A > small and T(A) > for A > is sufficiently 
smooth, we get that X(X) is sufficiently CS for A > small. The branch X(X) defined by 1)9.32(1 
is such that X(A)r(A) G QiD) for A > and it is the branch whose existence is stated in Theorem 

EH 

□ 

Proof of Corollary \5. 81 Using Proposition 15.71 we know that -'^(A) is sufficiently CS for A > o 
small. Therefore, |X(A)|^ for A > small is sufficiently CS. Using Taylor formula for e ^ |X(e^)| 
around e = up to the term of order e, and then taking e = y/X, we get |X(A)|^ = |X(0)|^ + O(e) = 
|X(0)|^ + 0{Vt) for A > smaU. Also, \X{X)\ is continuous for A > 0. 

The branch X-^(A) is constructed below. Let |Xo|T(0) = oq + 2/c7r, where oq S [0, 27r), G Z, 
/c > 0, and r(0) 



2lT 
l^bif I 



1. If \Xq\ 7^ kujfjif for all /c G Z, then oq ^ 0. 

Since A{T{0),0) = e^oT(o) ^ gX(o)T(o) ^ ^z(t,x) ^ ^^at 

X(0)r(0) = Z{t,0) = q{BCH{0,XoT{0))) = aoX^. Therefore, |X(0)| = > 0. Since X(A) 

is sufficiently CS, then |X(A)| > for A > small. 

Then Xf{X) = ^^^^^^^(A) is well-defined for A > smah. 

We check that X/(0) = Xq. Using X{0) = \X{0)\ X^, it yields Xf{0) = l^(o)ino)+2fc^ xi ^ 

ao+2fc-7r I V I Vl Y 

T(0) ^0 - \^0\ -^0 - -^0- 

Since |-^(A)| > and -^(A), ux are sufficiently CS, then X^{X) is sufficiently CS. 
Using the Taylor formula for e X^ (e^) around e = up to the term of order e and then 
taking e = A^ we get Xf{X) = X^(0) + 0{e)X^ + 0{e)Xi + 0{e)X2 = Xq + 0{X^)X^ + 
0{X^)Xi + 0{X^)X2 for A > smaU, thus we get the formula ((^ . 

We have only modulated rotating waves because |Xo|T(0) ^ 2kTT for any k €z Z implies 
|X-/'(A)| T(A) / 2kTr for any A; G Z and for A > small by the continuity of (A)| r(A). 



2. If \Xq\ = kuJbif for some k G Z, k ^ 0, then = 0. 
Let Xf{0) = Xq. 

If |X(A)| / 0, then we define X^X) = ^^^(^^^^(A) for A > small. 

If \X{X)\ = 0, then we define Xf{X) = X{X) + k \ujx\ QW, where Q{X) G so(3), |Q(A)| = 1 
for A > small. 
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If |Xo| = kojbif for some k e Z, k ^ 0, then we check that \Xf{X)\ = |X(A)| + |^ for all A > 

't(\)\X(X)\ 



small. For A > small such that \X{X)\ > 0, we have |X-^(A)| = \X{X)\ + 2A;7r^J^^^ 



\X{X)\ + j^. 

For A > small such that X{X) = 0, we have \Xf{X)\ = ^ \Q{X)\ = |gy. 
Hence, for A > smah, we get \X^ {X)\ = \X{X)\ + 

Since |X(A)|^ = \X{0)\^ + 0{VX) and |X(0)| = 0, we get \X{X)\ = 0{X^) and then |X-^(A)| = 
0(A3) + ^ + 0(A) = 0{X^) + \Xo\ for A > small, thus we get the formula 
Also, \X^ (A) I is continuous since X{X) and T(A) are continuous for A > small. 

By using Proposition (6) and (7)), we have e^'^^^^^^'^) = e^(^)'^(^) for A > small. □ 

Proof of Theorem \5.y[ 1. Suppose that |Xo| / kojuf for any k £ Z. 

We have 5^(t,A) = e-^^W*^(t,A) = e-^^W*e^(*'^) = e3(BCH(-x/{A)t,z(i,A))) ^j. ^ ^ p^^) 
and A > small. 

Since X-^ (A) and Z{t, A) are sufficiently CS for t G [0, oo) and A > small and the BCH 
formula in so(3) is smooth, it follows that Per-^{t,X) = q{BC H{—Xf {X)t, Z{t, X))) is suffi- 
ciently CS for t G [0, oo) and A > small. 

Since e-^'(0)*e^(*'0) = = h = e^CH{-xf (o)t,zm) ^ it follows that Per/(t,0) = 

g(SCif(-X/(0)t,Z(t,0))) = O3. Then, Perf{t,X) = \^Y{t,X) for all t G M and A > 
small. The periodicity of Y{., A) results from the following Remark: 

Remark 9.2. IfP{t,X) = e^^*^^\ P: MxM^ S0{2,) is a continuous T -periodic function in 
t for A > small and G: R x M — > so{3) is a continuous function such that G{t, A) = XH{t, X), 
then H{t, A) is a T-periodic function in t for A > small. 

Proof We have P{t + T, A) = P{t, A) =^ e^(*+^'^) = e^(*'^) =^ e^H{t+T) ^ ^XH(t) ^^^^ 
XH{t,X), XH{t + T, A) are in a neighborhood of O3 for A > small. Then, by using the 
fact that the exponential map exp is a local diffeomorphism at O3, it follows that XH{t, A) = 
XH{t + T, A) or G{t + T, A) = G{t, A) for A > small. □ 

Writing the Taylor formula for e Per^(t, e^) at e = up to order of e and taking e = \/A, 
we get Perf {t, X) = O3 + eYi{t, e) = y/XY{t, A)) for t £ [0, 00) and A > smaU. 

2. Suppose that \Xq\ = kuj^f for some A; G Z, /c 7^ 0. 

We have B{t,X) = e~^^^^* A{t,X) = e-^W*e^(*'^) = e9(^^^(-^W*'^(*.^))) for t G [0,oo) and 
A > small. Since X{X) and Z{t, X) are sufficiently CS" for t G [0, cxd) and A > small, and 
BCH formula in so(3) is smooth, it is clear that Per{t,X) = q{BCH{—X{X)t,Z{t,X))) is 
sufficiently C5 for t G [0, 00) and A > small. 

Since X(0) = O3, we get X(A) = 0{^/X)X^ + 0{^/X)Xl + 0{^/X)X2. Since ^(t,0) = e^o* = 
e^(*'°) and Per(t,0) = q{BCH{Os, Z{t,0))) = q{BCH{Os, Xot)), it yields that Per{t,0) = 
Xot{ mod 27r). We denote r(f) = t{ mod j^). Therefore, Per{t,X) = Xor{t) + VXHi{t,X) 
for t G [0, cx)) and A > small. Also, it is clear that B{t, A) = e^'>t+^/XH(t,x) _ g-^^^g ^(^^^ -g 

||^-periodic in t and P(0, A) = I3, it follows that e ^RI' = /g for A > smah. 

□ 
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Proof of Theorem mm Since ^{t,ux,X) = e^'-^^ B{t, X)^{q{t, X)) = e^^^^^ {t, X)^{q{t, X)), 
the conclusion of Theorem 15.111 result by applying Theorems 15. 7| 15.81 and 15.91 □ 



Proof of Theorem\7l\ Let us define Y{a,(3) = e"^i+^^2. 

Let us denote A = e^, X{e, fi) = X{e'^,fi) and T(e, jj) = T{e'^,fi) = y^^. Let r(0, /i) - 



We have to prove that there exist sufficiently smooth curves 

/.i = iJ-{e),a = a(e) and P = /?(e) such that /i(0) = 0, a(0) = and /3(0) = 

satisfying 

Y{a, (3) = e^(^'^)^(^'/^) or Y{a, /3)-ie^(^''^)^(^''^) = I3. (9.33) 
Let us define the function F: ^ 50(3) by 

F(a,/3,e,/x) = e-"^i-/3^2gX(e,M)f 

We will prove the existence of the smooth curves using the implicit function theorem for the function 
F at the point {a,(3,e,fi) = (0,0,0,0). 
We now prove that: 

1. F(0,0,0,0) =/3 
and 

2- [(-DF)(o,o,o,o)]xi,Xi,X2 rank 3. 



We have that F(0, 0,0,0) = e-o-^i-o-^2eX(o,o)T(o,o) ^ gX(o,o)T(o,o) ^ ^ \^,.f\ ^ ^^^^^ \Xo\ 
Recall that we have 

XG{q,X,fi) = xo{q,X,fi)X^ + xi{q,X,n)Xi +2 {q,X,fi)X2. 
Then Xg{0, 0, fi) = a{fi)X^ + b{ij)Xi + c{fi)X2, where 

a(/i) = xo{0, 0, fi) = \Xo\ + 0{n), b{fi) = xi(0, 0, fi) = 0{ijl) and c{^l) = ^2(0, 0, /x) = 0(/u). 

Since A(t,A,;u) = e^(^'^)^(^'^), then ^(t,0,/i) = e^^O-'^^^^O-^) = e^(0''^)^(°''^) . Also ^(t,0,/x) 
gXG{o,o,/i)r(o,/i)^ Using Proposition [Ol (2). we get 

_^ X(0,M)T(0,At)| ^ _ Ap^G(0,0,At)T(0,M)| 



„XG(o,o,o)r(o,o) ^ 

n=0 



00 r / -j^Njn 



^^^-^ (ad(XG(0,0)T(0,0)))" ^(XG(0,0,/i)r(0,^)) U=o 



(9.34) 
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We have r(0, 0) 



27r 



, Xg(0,0,0) =Xo and 



(XG(0,0,/x)r(0,/x)) 1^=0 = {XG{0,0,fx)) |^=oT(0,0) +Xg(0,0,0)(T(0,;u)) 1^=0 

\aiO)X^+b'iO)Xi + ciO)X2] 



\i^bif\ L 



a; 



27r 



(a'(0) - A: • (sgn(a;M/))^KM)'U=o)^o + b' {0)Xi + c'(0)X2 

(9.35) 



where we have used \Xq\ = kujbif- Therefore, H9.34() becomes 



1 T"^n 



n=0 



(-1)" / 27^\X^ 



(n + 1)! 



ad{- 



'■Ol ^1 
An 



^ . (a (0) - A:(u;o,/.)' U=o)^o + (0)^i + c (0)^2 



27r 



[(a'(0) - A:(a;o,M)'U=o)^o + ^'(0)Xi + c'(0)X2 



(a'(0) - fc(L^o,M)'U=o)^o + ^'(0)^1 + c'{0)X2 



^ (n + 1)! 

n=l 

^ [(a'(0) - k{u;oJ\f,=o)X^ + (6(0) + *)Xi + (c'(0) + *)X2 



(9.36) 







where we have used e l'^*''-''! = and where * denotes other terms. These term appear only in 
the coefficients of Xi and X2 since 



{ad{2\k\nXl)) 



27r 



27r 
4|A;|7r2 



\u!bif\ 

and for any integer n > 2 



(a'(0) - A:(c^o,m)'U=o)^o + ^-^O)^! + ^'(0)^2 
2\k\TTX^,{a\0) - A:(cuo,,.)'U=o)^o + ^'(0)^i + ^'(0)^2 
(-c'(0)Xi + 6'(0)X2), 



(9.37) 



(a'(0) - fc(u;o,^)'U=o)^o + &'(0)^i + c {0)X2 



2\k\7rXl (a(i(2|A:|7rXo^))" ' [(a'(0) - A:(a;o,M)'U=o)^o + &'(0)^i + c'(0)X2 



/ 27r 



(9.38) 



will contain only linear combinations of Xi and X2. 
Then with respect to the basis {Xq, Xi, X2}, 
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(^^)(0,0,0,0) 










-1 







-1 





V 



27r 
|l^M/ I 



(a'{0)-k{u;oJ\,=o) T^b' (0) + . rf^c'(0) + * 



I 



where * denotes other terms (Below we show that in fact (0) + * = and r^c' (0) + * = 0) 

^OO (-1)" i „^,2-K\Xa\ I 27T 



Letf/ = E-o^^arf(^^o^) 
ing Proposition IA.3I (1). we get 



(a'(0) - k{uJoJ\f.=o)X^ + b'{0)Xi + c'(0)X2 



Us- 



^ = E (ad(27r|fc|Xo^))" [(a'(0) - kM'\,=o)X', + 6'(0)Xi + c 



n=0 

27r 



(-ir 



(9.39) 



i2n\k\X^r f(o'(0) - A:(a;o,M)'U=o)^o + ^'(0)^1 + c'(0)X2 



If we apply the relation H9.22() for Z = 2TT\k\XQ, then taking into account that \Z\ = 2Tr\k\, we get 



E7^-^(2^i^i^or = ^3 + (xi)' 



n=0 



(n + 1)! 



(9.40) 



If we substitute (|On)) into (fOU]) . it follows that 
27r 



27r 



[/3 + [X^f] (a'(0) - A:(c^o,^)'U=o)^o + b'{0)Xi + c'(0)X2 



(a'(0) - k{LjQj\^=o)Xl + 6'(0)Xi + c'(0)X2 
+(a'(0) - fc(^o,M)'U=o)(^o')'^o + b'mx'.fX, + cmx'.fX, 
-Xi and {X^fX2 = -X2, we get U 



Since (Xo^)2xi = "o , {X^fX^ 



2-K 



(9.41) 



(a'(0) 



|fe|(u;o,^)'U=o)Xi. 

In order to have the rank of (^-^) (0,0,0,0) equal to 3, it is necessary and sufficient that 

a'(0) / A;(a;o,/,)'U=o- 

We note that depending on the function X^[t, X, n) we can establish if a(e) = and/or /3(e) = 0. 
Another proof of Theorem\TJ\We have that A{T{X, n), X, = e^MTM ^ where X(A,/i) is 
defined as in Proposition 15.71 

Let X(A,^) = ai(A,/u)X(J +6i(A,/i)Xi + ci(A,^)X2. Then, since X{X,^) is sufficiently CS, we get 

that ai, 61, ci are sufficiently CS and 6i(0,/i) = O{fj,),ci{0, fi) = 0{^). 

Since e^(0'0)^(°'°) = ^(r(0, 0), 0, 0) = I3, then X(0,0) = O3, that is ai(0,0) = 0. 

We want to find a sufficiently CS branch /i = /i(A) such that /i(0) = and ai(A,/i(A)) = 0, for 

A > small. If (ai)^(0, 0) 7^ 0, then by applying the implicit function theorem we get the existence 

of the required sufficiently CS branch. 
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Now we show that (ai)^j(0,0) 7^ is equivalent with a (0) ^ k{ujQ^^)' \^=q. 

We have that A{T{0, /x), 0, fi) = e^(0'^)'^(o,M) ^ gXG(o,o,M)T(o,M) ^j. |^| g^^^^^ ^-^^^ x{0, fj.)T{n, 0) 

Xg{0,0,ij,)T{h,0) + pff(5^n!^^G'(0'0"") '^^^^ ^ ^""^ small. 
Thus, [ai(0,^)Xi + 6i(0,//)Xi +ci(0,/i)X2] r(0,/i) = 

[a(^)Xi + + c(m)X2] r(0, fi) + [a(/")^o' + + c(/tx)X2] for > small. 

This implies that 

ai(0, //) = a{fi) + —=^^==l{n)uJo^f, • sgn{u;uf), (9.42) 

which for /i = gives = —\k\, where we have used the fact that |Xo| = kujhif > implies 
sgn{k) = sgn{u!hif), and thus, k ■ sgn{ujbif) = \k\. Since G Z is continuous, we get /(/u) = —\k\ 
for > small. 

Then by taking into account that sgn{ujf)if) = sgn{k), (|9.42p becomes 

ai(0, fi) = a(/i) - k "J^^ =u;o,f,. (9.43) 

By differentiation of 1)9. 43() . it follows that (ai)^(0, 0) = a'(0) — A;(ti;o,/i)'|/i=o / 0. 

□ 

Proof of Corollary \ 7. ^ The conclusions of Corollarv 17.21 result in the same way as in the proof 
of Theorems 15.71 and 15.91 except the scaling of the primary frequencies |X(A,^(A))| that results 
from Remark 3.5 in (53j. □ 



A BCH formula in so(3) 

It is known that 50(3) is diffeomorphic as a manifold to the real projective space MP'^ (see j42|). If 




c 



Y = —a c G •so(3), then we define ^ = b • Also, we define |y| 



Y 



A model 



a 



for the space MP^ is the set D = G | < vr, whith the antipodals points of the norm |y| = 
vr identified }. In fact, D is the quotient set -E/ ~, where ~ is the equivalence relation y ~ z iff z = 
~?/) \y\ = a-iid E = {~y £ M'^ | < vr}. The set D is considered with the quotient topology. 
Sometimes, instead of 'y' E M^, < vr, we use y G so{3), in which case we denote the equivalence 
class [y] = [Ij']. The projection map p: E ^ D , p{y) = [y] is smooth. There exists a unique 
smooth function d* : 50(3) D such that e'^*^^) = A. 

Definition A.l. For any X, Y e so{3), we define BCH{X,Y) = d*{e^e'^), where d* is defined 
above. Clearly, we have e^'^^(^''*^) = e^e^ and BCH{X,Y) G D, for any X,Ye so(3). 

Theorem A. 2 {BCH Formula in so(3), [4\ Ist l44] ) . The BCH formula in so(3) has the form 
BCH{X, Y) = [aX + pY + -f[X, Y]] for X,Y e so(3), (A.l) 

where 

a = k{X, Y)ha{X, Y),p = k{X, Y)hp{X, F), 7 = HX, Y)h^{X, Y), 
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and 



e = cos^cos^ -sin^sm^cos(Z(X,y)), 

. 1^1 \y\ 

ai = sm cos a = oie, 

oi = sm cos b = bie, 

. \X\ . \Y\ 

ci = sm sm , c = cie, 

2 2 



di = ^aj + 6f + 2ai6i cos(Z(X, y)) + ci{sm{Z{X, Y))^, 

d = di \e\ , 

where Z{X,Y) is the angle between the two vectors X and Y , 



ha{X,Y) 
hpiX,Y) 

hy{X,Y) = < 



cos^ lfX = 03, 



ifY^O^; 



m 

CosM ifY = 03, 



CI 



1^1 



1 



ifX^03,Y^03; 
tfX = 03,Y^03; 

ifY = 0^,X^Or, 
tfY = 03,X = 03, 



k{X,Y) 



arcsin(d) 
TT— arcsin{d) 



if (e^e^)^ 7^ /s, e^e^ has eigenvalues with positive real parts; 
if (e^e^)"^ 7^ Is, e^e^ has two eigenvalues with negative 
or zero real parts; 
^/(eV)2 = /3, e^e^^h; 
ife''e^ = h, 



where s = sgn{e) 



1, i/e>0; 
-1, i/e<0, 



for any {X,Y) G so(3) x so(3), such that (e"^e^)^ / h or 



I3. The functions o? , 0^ , 7^ are smooth on so{3) x so(3) and \a\, I7I are continuous 



e e 

on so{3) X so{3). Also, the function BCH is smooth from so(3) x so(3) into D. 
Proposition A. 3. 1. ^42] For any X, Y e so(3), 

adiX)Y = XY and 

X^Y for any integer n > 0. 



{ad{X)YY 

2. For any A G 5*0 (3) and any X G so(3), we have 



(A.2) 
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(a) 

\AXA~^\ = \X\; (A.3) 



AXA-^ = B if and only if aI = l3 . 
3. f2IF For any X £ so{3), we have: 

X2n ^ (-l)"-i|X|^(""^)x2/oranyn> 1, 



= /or any n>0. ^^'^^ 

The properties of the exponential map e: so(3) SO{3) are: 

Proposition A. 4. 1. The exponential map exp: so(3) — > SO{3) is surjective. 

2. J A' 61 \4<^ The exponential map exp is a smooth function on so{3) and its differential is given 
by: 

{d{exp))x{Y) = for any X,Yg so(3) 



n=0 

or 



(n + 1)! 



1 

((i(exi5))x(l^) = V -— — (ad(X))"(y)e^ for any X,Yg so(3) 

( we will use the first formula later ). Moreover, it is a local diffeomorphism near any 
X G so(3) if and only if the operator ad{X) has no eigenvalues of the form Inik with k ^ 0, 
that is if and only if \X\ ^ 2k'iT for k £ 1^, k ^ 0. 

^- m (e^^*'') = ^(t)e^(*) if and only if X{t) = Xg{t), where X G so{3) and g: R ^ M. is a 
function. 



4- 13 \44^ The Rodrigues' formula holds: 



2 \x\ 



V sin \X\ sin -^^ „ , , 

= I3 + ' + 2 ^X^ for any X £ so{3), 

\-^\ 

where we take the limit when X = O3. 

5. I42\l The exponential map exp maps any X € so{3), X ^ O3 to the right-handed rotation 
A G S0{3) with angle \X\ around X . 

6. 142^ e-^ = Is if and only if \X\ = 2k'K for some A; G Z. 

7. e^ = e^ if and only if either \X\ = 2/c7r = \Y\ orY=X + ^X for some keZ. 
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B Equivariant Center Manifold Reduction for Reaction-Diffusion 
Systems on rS^ 

We consider the reaction-diffusion system (|4.1jl . 

Theorem B.l ( [Ul |45j ) . Let uq £ Y" be a relative equilibrium that is not an equilibrium for \4-l\ ) 
at X = and such that the stabilizer of uq is S„q = I3. Let L be the linearization of the right-hand 
side of with respect to the rotating wave uq, 0) = e^°^Uo at X = in the co-rotating frame, 
that is 

L = DAs + DuF{uo,0) - Xo. 

Suppose that: 

1. (y{L) n {z G C I Re [z) > 0} is a spectral set with spectral projection P*, and dim{R{P^:)) < 00; 

2. the semigroup e^* satisfies |e^*|R(i-p,) | < Ce"^"* for some /3o > and C > 0. 

Then there exists a sufficiently smooth parameter- dependent center manifold M^"(A) for the relative 

equilibrium uq. Let Kk be the orthogonal complement of Tuq{S0{3)uq) in = R{P^:). 

Then the center manifold M^"(A) is diffeomorphic to S0{3) x Kk for \X\ small. Furthermore, there 

exist sufficiently smooth functions Xq : T4 x M — > so(3) and X^ : x M ^ Kk such that any solution 

of 

A = AXG{q,X), /T.^^ 
q = Xr,{q,X), 

on S0{3) X Kk corresponds to a solution of the reaction- diffusion system J^. j| ) on M^"(A) under 
the diffeomorphic identification {A,q) A^{q), with ^ a local chart from V to M, for \X\ small. 
Also, ^^(0, 0) = Xq, X]\f{0,0) = and a{Dy,Xpf(0,0)) = (t((5*L|v;), where is the projection 
onto along Tuq{SO{3)uo). 



Remark B.2. Theorems \BA\ is valid for more than one parameter. 

C Some Computations 

We present tlie proof of ()8.1|) : 

If we parameterize 5*0(3) by Euler angles, that is j4 = Rz{'4')Rx{S)Rz{4')i where (j) ^ [Oj^tt), 
6* G [0, vr], ^ G [0, 27r), and we substitute this into the equation A = AXciq, A), we get 

= R,{ij)Rx{e)R,{4>) [F^{q, X)Lx + Fy{q, X)Ly + F'^iq, A)L,] (C.l) 



or 



i^R,{-(t))Rx{-d)R,{-^) (^(V^)) Rz{(p)RMRz{<t>) 

+ eR.{-4>)R.{-e) (^(^)) R.{<t>) + ^Rz{-4>) (^i^)^ ^^'^^ 

= F^{q, X)L., + Fy{q, X)Ly + F'{q, X)L,. 
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By computation, it follows that 



R,{-mx{-e)R^{-^) (ff (V')) RM)Rx{G)Rz{<t^) 

If we substitute (l(T3|) into (l(T2l) . we get 



= (— sin + cos (/iL^., 

= cos + cos sin 0Lj^ + sine/) sin 0La;. 

(C.3) 



ijj (cos OLz + cos 4> sin 6Ly + sin </> sin 6*^^ ) + ( ( — sin (j))Ly + cos (/jL^,. ) + (j)Lz 

= F^q, X)L, + Fy{q, \)Ly + A)L,. (C.4) 

If we identify the corresponding coefficients of L^, Ly and from both sides of the equation HC.4|) . 
we get the system 

'ijjcose + (j) = F^q,X), 
tpcos(f)sme - 6sm(p = Fy{q,X), (C.5) 
tp sin (j) sin 9 + 9 cos cp = F^ (9, A) . 

Then, by multiplying the second equation in ()('.5|) by — sin cj) and the third equation in ()(L5|) by 
cos (j) and adding them, we get = —F^^q, A) sin (j) + F^{q, A) cos (/>. Similarly, by multiplying the 
second equation in ()('.5|) by cos (p and the third equation in ()('.5|) by sin (p and adding them, we 
get Tp = [F'y{q, A) cos (f> + F^{q, A) sin 0] . If we substitute ip given by the above relation into the 
first equation of the system HC.5|) , we get (p = F^ {q, A) — cot 9 A) cos (p + F^{q, A) sin (p] . 
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